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ABSTRACT 

We present a method which uses colour-colour cuts on SDSS photometry to select 
white dwarfs with hydrogen rich (DA) atmospheres without the recourse to spec- 
troscopy. This method results in a sample of DA white dwarfs that is 95% complete 
at an efficiency of returning a true DA white dwarf of 62%. The approach was applied 
to SDSS Data Release 7 for objects with and without SDSS spectroscopy. This lead 
to 4636 spec troscopicially con fi rmed DA white dwarfs with g ^ 19; a ~ 70% increase 
compared to Eisenstei n et alj 's I2Q06I sample. Including the photometric-only objects, 
we estimate a factor of 3 increase in DA white dwarfs. We find that the SDSS spectro- 
scopic follow-up is 44% complete for DA white dwarfs with T c ff > 8000K. We further 
cross-correlated the SDSS sample with Data Release 8 of the UKIDSS Large Area 
Survey. The spectral energy distributions of both subsets, with and without SDSS 
spectroscopy, were fitted with white dwarf models to determine the fraction of DA 
white dwarfs with low-mass stellar companions or dusty debris discs via the detection 
of excess near-infrared emission. From the spectroscopic sample we find that 2.0% of 
white dwarfs have an excess consistent with a brown dwarf type companion, with a 
firm lower limit of 0.8%. From the white dwarfs with photometry only, we find that 
1.8% are candidates for having brown dwarf companions. Similarly, both samples show 
that ~ 1% of white dwarfs are candidates for having a dusty debris disc. 

Key words: Stars: individual: white dwarfs - infrared excess 



1 INTRODUCTION 

White dwarfs are the most common stellar remnants in 
the Galaxy, having descended from main sequence stars 
with 0.8M Q <M<8M©. Of these, hydrogen rich (DA) 
white dwarfs make up the vast majority. Large sam- 
ples of white dwarfs are particularly useful in many 
studies, for example constraining the luminosity func- 
tion, which in turn can be used to determine the ages 



can probe the companion mass function over a wide range of 
(initial) host star masses. Currently, only four white dwarfs 
are confi rmed to have (non-interacting) substellar com- 
panio n s dBecklin fc Zuckermanl 19881: Farihi fc Christopher! 



of many Galactic populations (e.g 



lOswalt et all 1 19961 : |Pe Gennaro et all 



Winget et all 1 19871 : 
20081 ). The low lu- 



| 2004 iMaxted et al.l 120061 : ISteele et alj 120091 . but see 
iLuhman et al.ll201ll for a very low-mass candidate) , and the 
fraction of white dwarfs with brown dwarf companions ap- 
pears to be consistent with the low number fou nd around 
FGK stars (|Farihi et al.ll2005l : iHoard et al.ll2007l h 



minosity of white dwarfs also makes them ideal targets 
for searches of low-mass c ompanions, such as pioneered by 
IProbst fc Oconnelll l|l982h . The companion mass distribu- 
tion is thought to drop near the low- mass end (~ O.IMq), 
and the fraction of FGK stars with substellar compan- 
ions estimated from radial velocity surveys is < 1% (e.g. 
Marcv fc Butlerll2000l: iGrether fc Lineweave3l2006h , though 
Metchev fc HillenbrandT (l2009 ) suggest that substellar com- 
panions are more frequent at larger orbital separations. Be- 
cause white dwarfs are the progeny of main-sequence stars 
with masses of up to 8 Mq , studies of white dwarf binaries 



While white dwarfs are a lso excellent candidates for 
searc hes of Jovian planets (e.g. llgnacdfeOQll ; iBurleigh et al.l 
120021) ■ no planet has y et been unambiguous l y discovered 
around a white dwarf (IBurleigh et all 120081: [Hogan et all 
120091 : iMullallv et alJ[200gf T However, ~ 15 white dwarfs arc 
known to exhibit infrared flux excesses that can not be ex- 
plained by the p resence of low-ma ss (planetary or substel- 
lar) companions (|Farihi et alj 2009). The first such infrared 
excess was found around G 29-38 a nd tentatively associ- 
ated with a brown dwarf companion l|Zuckerman fc Becklinl 
1987). Using the fac t that G 29-38 is a pulsating white 
dwarf, iGraham et"aH l|l990h demonstrated that the obser- 
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vations are best explained with the presence of circum- 
stellar dust heated by the white dwarf, and proposed the 
disruption of a larg e asteroid as the origin of the dust. 
iKoester et al.l |l997h detected metals in the photosphere of 
G 29-38, demonstrating that the white dwarf is accreting 
circumstellar material. The scenario of the t i dal di sruption 
of rocky asteroids has been refined by ljural (|2003l ) , and is 
now the generally accepted interpretation of the growing 
number of metal-p olluted white dwarfs with infrared excess. 
iFarihi et al.l |2009l ) estimate that 1 — 3% of all white dwarfs 
with cooling ages below 0.5 Gyr should exhibit infrared ex- 
cess related to the tidal disruption of asteroids. The location 
of the debris material within the white dwarf tidal disruption 
ra dius for a rocky astero id has ultimately been confirmed 
bv lGansicke et all ((2006b) who detected Doppler-broadened 
Ca II emission lines from a gaseous component of the circum- 
stellar disc in SDSSJ122859.93+1 04032.9, a metal-pollut ed 
white dwarf with infrared excess (|Brinkworth et aflljooj ). 

In this paper, we use the Sloan Digital Sky Survey 
(SDSS) to first identify a method to select hydrogen-rich 
(DA) white dwarfs using photometric criteria, and then de- 
termine those with infrared flux excess in the UKIRT In- 
frared Deep Sky Survey (UKIDSS). Our method is sensitive 
to unresolved M and L-type companions and to warm de- 
bris discs. The long-term goal of this project is to increase 
the number of such objects available for detailed follow-up 
studies, and to provide a more detailed understanding of the 
frequency of both substellar companions to and planetary 
debris discs around white dwarfs. 

The structure of the paper is as follows: In Sect.[2]we 
discuss the three surveys that are utilised in this paper; 
the Sloan Digital Sky Survey, the UKIRT Infrared Deep 
Sky Survey and the Wide-field Infrared Survey Explorer 
(WISE). The method for selecting DA white dwarfs from 
SDSS DR7 is outlined in Sect .[3] along with a discussion 
of the spectroscopic completeness of SDSS with respect to 
DA white dwarfs. The procedure for cross-matching ob- 
jects from SDSS with UKIDSS is demonstrated in Sect. [4] 
Identifying which white dwarfs show an IR excess and the 
method for fitting these excesses is given in Sect. [5] Further 
to this, Sect. 15.41 discusses comparisons between the spec- 
troscopic and photometric fitting methods, including how 
contaminants were dealt with. A summary of the overall 
numbers of objects at each stage of the process is given in 
Sect. 15.51 Some notes on interesting objects discovered in 
this work are given in Sect .[6] A comparison between the 
work he rein and the white dwarf-main sequence star cata- 
logue of iRebassa-Mansergas et all l|2010l ) is made in Sect. [3 
Objects that are detected in WISE are discussed in Sect. [8] 
The implications of the DA white dwarf selection and the 
IR excess selection methods are discussed in Sect .[9] Finally, 
the conclusions of this work are drawn together in Sect . 1101 




0° ' 90° 18*0"° 270' 360° 

R.A. [degrees] 



Figure 1. Coverage of the SDSS DR7 (grey) and UKIDSS DR8 
LAS (black) in equatorial coordinates. 

2.1 SDSS 

We have made use of the SDSS Data Release seven (DR7, 
lAbazaiian et alj 120091 ). which represents the final DR of 
the SDS S II project, inclu ding the low-latitude extension 
SEGUE (|Yannv et al.ll2009l ). SDSS DR7 provides ugriz pho- 
tometry for 357 million objects spanning a magnitude range 
~ 15 — 22 and covering 11500 deg 2 , approximately one- 
quarter of the celestial sphere, as well as follow-up low- 
resolution [R ~ 1800, 3800 - 9200 A) spectroscopy for 1.44 
million galaxies, quasars, and stars. We limit our search 
to point-like sources, and to apparent magnitudes g ^ 19, 
which leaves ~ 24 million unique photometric and 87,000 
unique spectroscopic objects. An additional important re- 
source within DR7 are proper motions computed from the 
USNO-B and SDSS positions. 

2.2 UKIDSS 

UKIDSS is a set of five near infrared surveys being un- 
dertaken with the Wide Field Camera (WFCAM) instru- 
ment on the United Kingdom Infrared Telescope (UKIRT) in 
Hawaii. One of the five sub-surveys, the Large Area Survey 
(LAS), aims to be the IR counterpart to the SDSS. Here we 
match SDSS hydrogen-dominated (DA) white dwarfs with 
UKIDSS LAS as opposed to 2MASS used in similar stud- 
ies (e.g. iHoard et al.l 120071 ) because a large majority of the 
SDSS white dwarfs are too faint in the IR for 2MASS to 
pick up. UKIDSS LAS will eventually provide imaging over 
4028 deg 2 in four broad band colours, Y, J, H, and K, with 
limiting (Vega) magnitudes of 20.2, 19.6, 18.8 and 18.2, re- 
spectively, which adds a significant increase in depth over 
2MA SS. Here, we made use of UKIDSS DR8 (see lDve et all 
120061 ). The overlap between SDSS DR7 and UKIDSS/LAS 
DR8 is illustrated in Fig.[T] and amounts to ~ 2700 square 
degrees. 



2 LARGE AREA SURVEYS 



2.3 WISE 



The Sloan Digital S ky Survey (SDSS, lYork et al l l2000l : 
IStoughton et al] 120021) and the UKIRT Infrared Deep Sky 
Survey (UKIDSsT lLawrence et all l2007h are currently the 
deepest large-area optical and infrared surveys that are pub- 
licly available, and we summarise below details of both sur- 
veys relevant for our work. 



The Wide-field Infrared Survey Explorer (WISE) is a NASA 
Medium-class Explorer mission designed to survey the en- 
tire sky in four infra red wavelengths, 3.4, 4.6, 12, and 22/im 
I Wright et al] l201fj| ). WISE consists of a 40 cm telescope 
that images all four bands simultaneously every 11 s. It cov- 
ers nearly every part of the sky a minimum of eight times, 
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ensuring high source reliability, with more coverage at the 
ecliptic poles. Astrometr ic errors are less t han 0.5 arcsec 
with respect to 2MASS l| Wright et al.ll20ld ). The prelimi- 
nary estimated S /N = 5 point source sensitivity on the eclip- 
tic is 0.08, 0.1, 0.8, and 5 mjy in the four ban ds (assuming 
eight exposures per band; IWright et al ] |2010h . Sensitivity 
improves away from the ecliptic due to denser coverage and 
lower zodiacal background. We took advantage of the pre- 
liminary data release (PDR). The SDSS white dwarfs that 
are the subject of this work are largely too faint to be seen 
in the 12, and 22^m bands and therefore we generally only 
discuss the 3.4 and 4.6^im fluxes. 



3 SELECTING DA WHITE DWARFS IN SDSS 
DR7 

The latest catalogue of spectrosco picially identified SDS S 
white dwarfs was based on DR4 (jEisenstein et all [2006). 
which comprised roughly half of the sky coverage of DR7. 
Here, we exploit the much larger footprint of SDSS DR7, 
and also extend the white dwarf sample to photometric 
objects without follow-up spectroscopy. We restricted our 
ambitions to DA white dwarfs for a number of reasons. 
Firstly, the vast maj ority of all known wh ite dwarfs be- 
long to the DA class i|McCook fc Sionlll999l) . Secondly, de- 
termining the atmospheric parameters of DA white dwarfs, 
T c ff and log(gr), fr om fitting atmospher e models to ei- 
ther spectroscopy dBergeron et al.l 1 19921 ) or photometry 
l|Koester et al.lll979l ) is a well-established and robust pro- 
cedure. This is essential for the purpose of identifying white 
dwarfs with infrared excess, as we need to accurately ex- 
trapolate the white dwarf flux to the J, H, and K bands. 
Thirdly, optical spectra of DA white dwarfs are characterised 
by strong Balmer absorption lines on a blue continuum, and 
the strong dependence of the Balmer line equivalent widths 
results in DA white dwarfs occupying a distinct region in 
colour space. 

We have developed a two-pronged approach to iden- 
tify as many DA white dwarfs with spectroscopy within 
DR7, and subsequently to select white dwarf candidates 
which have ugriz photometry but were not spectroscopi- 
cally followed-up by SDSS. 

As a start, we retrieved the DR7 spectra and ugriz 
photom etry for all whit e dwa rfs with g ^ 19 and classi- 
fied by lEisenstein et al.1 (|2006h as DA or DA_auto, corre- 
sponding to visually confirmed and automatically classified 
hydrogen-dominated white dwarfs, respectively. This totals 
2889 unique objects, 938 being classified as DA and 1951 
DA_auto. All spectra were visually inspected to corrobo- 
rate their DA classification, and we found 99.4% agreement 
with t he classification for white dwarfs by lEisenstein et al.l 
(2006). The 0.6% disagreement primarily comes about from 
non-DA white d warfs that were classified as DA_auto by 
lEisenstein et alj 's classification routine. The sample of spec- 
troscopically confirmed DA white dwarfs was then used to 
trace the locus of DA white dwarfs in the (u — g,g — r), 
(g — r,r — i), and (r — — z) colour-colour planes. The pop- 
ulation of DA white dwarf follows a boomerang-shape in 
(u — g,g — r) colours, where it is clearly separated from the 
main sequence, but intersects the quasar population (Fig. [2}. 

This relatively complex shape was approximated by 



the intersection of two 5th-order polynomials (Table[T}. In 
(g — r, i — r) the DA white dwarfs lie along a relatively narrow 
band, overlapping with the blue end of the main sequence 
and, to some extent with quasars, which we approximated by 
the combination of a 2nd and 3rd order polynomial. Finally, 
in (i — i, i — z) , the DA white dwarfs are again located at the 
blue end of the main sequence, but display a relatively large 
spread in i — z. Visual inspection of the SDSS spectra of the 
outliers confirms them as DA white dwarfs, though the ma- 
jority of them are near the faint end of the lEisenstein et al.l 
sample. We decided to include those outliers in our colour- 
cut, and hence approximated the DA locus by the intersec- 
tion of three linear relations in (r — i, i — z). 

We applied a magnitude cut of g ^ 19 as a (conserva- 
tive) measure, so that each of the SDSS white dwarfs within 
the UKIDSS footprint would have a significant A"-band de- 
tection. Finally, we applied a set of recommended data qual- 
ity flags to the SDSS photometry to minimise contamination 
by instrumental artifacts and blended stars. Applying the 
constraint set summarised in Table[T]to DR7 resulted in the 
selection of 7444 unique spectroscopic objects, which were 
then visually classified as DA white dwarfs, other (non-DA) 
white dwarfs, quasars, narrow line hot stars (see below), and 
other objects (Tabled • 

3.1 Narrow Line Hot Stars 

The optical spectra of ultra-low mass white dwarfs and very 
cool DAs, particularly those at low signal-to-noise ratio, can 
look rather similar to early type main sequence stars, sub- 
dwarfs, extreme horizontal branch stars or very metal-poor 
halo stars, which we all lump together as contaminants with 
the designation "narrow line hot stars" (NLHS). These are 
particularly dominant at the bright end of the white dwarf 
banana (g ^ 16, see Table[2]). While ther e is noticeable inter- 
est in ultra-low mass white dwarfs ( e.g. iLiebert et al. I l2004l ; 
iKilic et al.ll2007l ; iMarsh et aill2010l ). they represent a tiny 
fraction of all DAs. Given that our aim is to study the bulk 
population of DA white dwarfs, we make no attempt to ac- 
curately classify ultra-low mass white dwarfs. In addition, 
our colour selection is only complete to effective tempera- 
tures above ~ 8000K. Attempting to include much cooler 
white dwarfs would result in significant contamination from 
NLHS. 



3.2 Overall Completeness and Efficiency 

Completeness and efficiency were the key parameters in de- 
signing our selection algorith m (Table[p, wher e we define 
completeness as the fraction of lEisenstein et ail 's DA white 
dwarfs recovered by our constraints, and efficiency as the 
ratio of spectroscopically confirmed DA white dwarfs to all 
objects in our colour-magnitude selection. In the case of the 
spectroscopic DA sample, one may argue that completeness 
is more important than efficiency because contaminants can 
be removed through visual spectral classification of all candi- 
date objects. On the other hand, the photometric-only DA 
sample requires a high level of efficiency to minimise the 
number of contaminants. We optimised the colour bound- 
aries to maximise both completeness and efficiency, and the 
constraint set in Table[T] results in a completeness of 95.4% 
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Table 1. Colour selection for finding DA white dwarfs in ugriz space. Objects were selected to be primary objects and point sources. 
Flags are shown in Hexagesimal notation. These are the standard SDSS "good photometry" flags, as documented on SDSS CAS JOBS 
jETfe Thakad 120081) , and were chosen so that the object was detected in BINNED1, and did not have any of the following: EDGE, 
NOPROFILE, PEAKCENTER, NOTCHECKED, PSF.FLUXJNTERP, SATURATED, or BAD.COUNTS.ERROR. 



Colour 


constraint 






(« ~ 9) 




-20.653 x (g - rf + 10.816 X (g - rf + 15.718 x (g 


- rf - 1.294 x (g - rf - 0.084 x (g - r) + 0.300 


(« - 9) 




-24.384 x (g - rf - 19.000 x (g - rf + 3.497 x (g - 


rf + 1.193 x {g - r) 2 + 0.083 x (g - r) + 0.610 


(9 ~ r) 




-0.693 X (r - if + 0.947 X (r - i) + 0.192 




(9 - r) 




-1.320 X (r - if + 2.173 X (r - if + 2.452 X (r - i) 


- 0.070 


(r — i) 




-0.560 




(r — i) 




0.176 x (i - 2) + 0.127 




(r — i) 




-0.754 x (i- 2) + 0.110 




9 




19 







j = 


flags & 0x10000000 









flags & 0x8100000c00ac 






Figure 2. Colour-colour plots illustrating the location of the SDSS spectroscopic objects. DA white dwarfs, non-DA white dwarfs, 
quasars and main-sequence stars are shown as blue, magenta, green and grey dots respectively. The polynomial colour cuts from 
Table[T]are overlaid as red lines. 



and an efficiency of 62.3% (Tabled)). From SDSS DR7, a 
total 4636 unique spectroscopic DA white dwarfs are con- 
tained within the colour-magnitude cuts. This represents a 
70% incr ease in spectrosc o pic D A white dwarfs with (7 ^ 19 
over the Eis enstein et al.l l|2006r ) sample. The photometric- 
only DA white dwarf candidates sample similarly contains 
9341 objects. Assuming 62.3% efficiency of the selection, 
~ 6000 extra white dwarfs are cont ained within the s ample 
and therefore the total increase over lEisenstein et al l (|2006l ) 
is approximately a factor or three. The efficiency is how- 
ever only a lower limit for the photometric sample because 
the SDSS is by design almost spectroscopicially complete for 
QSOs, one of our main contaminants. Therefore, QSOs are 
only a minor contaminant in the photometric-only sample. 

The spectral classification and completeness are given 
as a function of p-magnitude in Tabled The completeness 
drops slightly towards larger apparent magnitude because of 
the larger scatter in the colour-colour diagrams. The frac- 
tion of NLHS contaminants is largest at the bright end of 
our sample, which is a natural consequence of the much 
larger intrinsic brightness of subdwarfs and early-type main- 
sequence stars, and the fraction of quasar contaminants 
increases towards larger apparent magnitudes. It is possi- 
ble to eliminate a fair fraction of the contaminants in the 



photometric-only sample by using additional information 
such as proper motions and infrared colours (see Sect. 15.41) . 

3.3 Completeness of SDSS spectroscopy for DA 
white dwarfs 

The sample produced here provides an excellent opportu- 
nity to investigate the spectroscopic completeness of SDSS 
for DA white dwarfs. We used the cuts in Table[T]to select 
both the spectroscopic and photometric objects within DR7. 
From these two sets of data, we then calculate the spectro- 
scopic completeness within the (it — g,g — r) colour-colour 
plane (Fig.[3]). The upper middle and right panels display 
the density of spectroscopically confirmed DA white dwarfs 
(middle) and contaminants (such as NLHS and quasars, 
right hand side panel). For comparison, the upper left 
hand panel shows th e DA white dwarf cooling tracks from 
iHolberg fc Bergeron! (|20Q6TFl The efficiency of our colour 
cuts is obtained for each bin within the (u — g,g — r) plane 
as the ratio of the number of the DA white dwarfs to the 
total number of objects in the bin (lower centre panel). This 

1 See http:/ /www. astro. umontreal.ca/~bergeron/CoolingModels 
for an updated grid. 
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Table 2. Efficiency and completeness of the polynomial colour-colour cuts (Table[TJ as a function of SDSS g magnitude. The total number 
of spectroscopic objects that matched our selection and were visually classified was 7444. The classification 'Other white dwarfs' contains 
white dwarfs of type DAB, DAO, DB, DC, DQ, DZ, magnetic white dwarfs and white dwarf+MS binaries. The classification 'Other' 
includes CVs, galaxies, peculiar objects and unidentified spectra. The final two right hand columns show the number of photometric-only 
objects in each g magnitude bin (Nx t) an d the number of predicted DA white dwarfs, calcul ated by multiplying this by the efficiency 
(Nda)- The bottom line of the table gives the total number of DA white dwarfs with (g ^ 19) in Eiscnstcin ct al. 1 2 0061 ). and the number 
of them included in our colour cut. 





All WD 


DA WD NLHS QSO 




Other WD 


Other 


Efficiency 


Completeness 


Photometric-only 


9 


Candidates 


N % N % N 


% 


N % 


N % 






N To t 


N DA 


< 16 


258 


79 30.6 161 62.4 


0.0 


13 5.0 


5 1.9 


30.6% 


95.8% 


893 


273 


16-17 


581 


326 56.1 185 31.8 22 


3.8 


44 7.6 


4 0.7 


56.1% 


96.7% 


980 


550 


17-18 


1719 


1092 63.5 225 13.1 230 


13.4 


165 9.6 


7 0.4 


63.5% 


96.5% 


2278 


1447 


18-19 


4886 


3139 64.2 269 5.5 1028 


21.0 


439 9.0 


11 0.2 


64.2% 


95.0% 


5190 


3332 


Total 


7444 


4636 62.3 840 11.3 1280 


17.2 


661 8.9 


27 0.4 


62.3% 


95.5% 


9341 


5819 




Total DA WD 


Photometrically selected DA 


s 








Completeness 






Eisenstein 


2889 


2757 










95.4% 







clearly displays a reduced efficiency of selecting both the 
hottest and coldest white dwarf because of the increased 
numbers of contaminants. Our selection method however 
has an extremely high efficiency when selecting white dwarfs 
with temperatures between ~ 10, 000 — 20, 000K. The num- 
ber of DA white dwarfs without SDSS spectroscopy is pre- 
dicted by scaling the number of photometric-only objects 
with the efficiency (resulting in the lower left panel). Fi- 
nally, the DA white dwarf spectroscopic completeness was 
calculated as the ratio of spectroscopically confirmed DA 
white dwarfs to the total number of DA white dwarfs, with 
and without spectroscopy (lower right panel). The over- 
all spectroscopic completeness is 44.3% down to g — 19. 
As mentioned in Sect . 13.11 this analysis is limited to white 
dwarfs with T off > 8000K. The preference of SDSS spec- 
troscopy to target quasars is clearly seen in the lower right 
corner of the spectroscopic completeness colour-colour di- 
agram, where the vast majority of quasars lie. This corre- 
sponds to a very high spectroscopic completeness for cool 
DA white dwarfs. In contrast, the spectroscopic complete- 
ness for white dwarfs with T c ff > 12000K is significantly 
lower. 



4 CROSS-MATCHING WITH UKIDSS 

All (spectroscopic and photometric-only) objects from the 
DA selection in SDSS DR7 were matched with the UKIDSS 
database using the CrossID function. To decide upon a 
matching radius, a sample of 5000 randomly selected spec- 
troscopically confirmed DA white dwarfs were matched to 
the UKIDSS database with a 60" search radius r. The dis- 
tribution of the distance between the SDSS objects and the 
UKIDSS matches is shown in Fig. 2] The number of all pos- 
sible matches within 60" (black histogram) grows approx- 
imately as r 2 , as would be expect for chance coincidence, 
whereas true matches are primarily within r < 3". Select- 
ing only the closest match (blue histogram), the majority 
of these random mismatches are removed. The blue and 
black distributions agree well at small distances (r < 2.5"), 
indicating that crowding is not a major problem. 2.5" is 
much larger than the quoted astrometric accuracies of both 



Table 3. The number of all SDSS DR7 objects satisfying our con- 
straint set (TablcJTJ, and of various subsets with different UKIDSS 
bands. 



Detections 


Spectroscopic 


Photometric 


in band 


Objects 


Objects 


Total SDSS 


7444 


9341 


Any UKIDSS 


1990 
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1815 


1614 
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1549 
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1503 


1281 
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1108 


840 


H & K 


1075 


809 


J,H & K 


979 


720 



SDSS and UKIDSS (of order a few tenths of an arc second), 
but the large proper motions of the white dwarfs and the 
potentially large (~few years) time interval between both 
surveys can lead to positional shifts up to a few arc sec- 
onds. We adopt r — 2.5" for the final cross-matching of 
our spectroscopic and photometric-only SDSS DA samples 
with UKIDSS, which limits the number of spurious matches, 
and will exclude only a handful of (halo) objects with ex- 
tremely high proper motions. Any remaining positional mis- 
matches are flagged in the examination of the SDSS and 
UKIDSS images carried out later. We restrict our analysis 
to unresolved systems, as the physical association of spa- 
tially resolved companions to white dwarfs will be difficult 
to demonstrate with the available data. Consequently, ob- 
jects that are flagged as partially resolved in the UKIDSS 
images are removed form the sample. 

A total of 1990 of the SDSS objects with spectra 
were found to have at least one measured magnitude in 
the UKIDSS database. 1275 of these are spectroscopi- 
cally confirmed DA white dwarfs. Similarly, 1771 of the 
photometric-only objects had at least one match in the 
UKIDSS database. 
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Figure 3. The spectroscopic completeness of DA white dwarfs in SDSS DR7 within the (u — g,g — r) colour-colour plane. In the 
to p left panel, the colo u r sele ction from Table[TJ is shown as a red dashed line, overlaid with the DA white dwarf cooling tracks 
of Holbcrg & Bergeron (2006). From the bottom up, these curves represent log(g) = 7.0 — 9.5 in steps of 0.5. The number of 
spectroscopically confirmed DA white dwarfs and contaminants within our colour selection are shown in the upper middle and 
right panels, respectively. Two distinct regions of high contamination are visible, with NLHS and quasars being concentrated at 
the bluest and reddest colours of the DA "banana" , respectively. The number of DA white dwarfs without SDSS spectroscopy is 
shown in the bottom-left panel, calculated as the number of photometric-only objects weighted by the colour-dependent efficiency 
of our selection algorithm (lower middle panel). Finally, the spectroscopic completeness of SDSS for DA white dwarfs, i.e. the 
ratio of spectroscopic DAs to the total number of DAs, is shown in the bottom right panel. Cool white dwarfs have a very high 
completeness thanks to their colour- proximity to ultraviolet-excess quasars, which were intensively targeted by SDSS. 



5 IDENTIFYING INFRARED EXCESS 
OBJECTS 

5.1 DA White Dwarf Fitting 

A grid of synthetic DA white dwarf spectra was cal- 
culated with t he model atmosphere code described by 
Kocster and using the latest line profiles of 

Tremblav fc Bergeron! <|2009h . These cover T eff = 6,000 - 
100, 000K in 131 steps nearly equidistant in log(T e ff), and 
l°g(3) = 5.0 — 9.5 in steps of 0.25dex. 



5.1.1 Fitting the SDSS spectroscopy 

We fitted the SDSS spectra of all DA white dwarfs found 
within ou r colour cuts (T ablcPJ following the method de- 
scribed in lRebassa-Mansergas et al.l (|2007l ). A x 2 minimisa- 
tion is used to find a best fit from our grid of DA white dwarf 
model spe ctra, providing log(ff) and T e ff. Using the cooling 
models of iHolberg k, Bergeron! (|2006l 'l . M wd , R w d and the 
distance d can be calculated for each object. 



5.1.2 Fitting the SDSS photometry 

We also fitted all photometric objects found with our colour 
cuts, including all objects that do have SDSS spectroscopy. 
We also fitted objects known not to be DA white dwarfs, to 
allow us to investigate the properties of the contaminants 
among the photometric-only DA candidates. 

Photometric objects were fitted by comparing the SDSS 
u, g, r and i magnitudes to the white dwarf model grid, again 
based upon a smallest % 2 ■ Model magnitudes were calculated 
for each of the log(g) and T e fi values by folding the theoret- 
ical spectra through the SDSS and UKIDSS ugrizYJHK 
filter curves. Magnitudes redder than i were not included 
since we are searching for composite systems. This ensures 
we can recognise objects with an excess already showing 
in 2, such as white dwarf plus M-dwarf binaries. This does 
not significantly affect the spectroscopic fitting method be- 
cause we fit line profiles. This method is most sensitive at 
shorter wavelengths, where the companion does not signif- 
icant contribute. For the majority of photometric DA can- 
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Figure 4. Spatial offsets of the SDSS and UKIDSS positions for 
a sample of 5000 spectroscopically confirmed DA white dwarfs 
randomly selected from our constraint set (Table [2}. The blue 
histogram shows the distance to the closest neighbour in UKIDSS, 
the black histogram plots the distances to all possible matches. 
The bins have a width of 0.1". 




1.0 1.5 
Wavelength [pm] 



Figure 5. The SED of SDSS J1218+0042, an example of a spec- 
troscopically confirmed DA white dwarf with a possible IR flux 
excess (classified as DAire: in Table 0, SDSS J1218+0042. The 
best-fit model to the ugri photometry is shown (T c g = 10000K, 
log(g) = 7.75), which results in an H-band excess just below 3a. 
Adopting the slightly higher temperature (T c g = 11173K) from 
the analysis of the SDSS spectrum boosts the excess to just above 
3c, flagging the object as an IR excess candidate. However, given 
that Y and J-band magnitudes fall significantly below the model, 
further IR data are necessary to confirm or refute the IR excess 
of this white dwarf. 



didates, the 4-band photometry did not provide sufficient 
constraints to accurately determine the surface gravity. For 
objects with calculated effective temperatures in the range 
9000 — 20000K, where the spectral line widths are narrow 
and therefore do not significantly affect the SED, we adopted 
a canonical value of log(<;) = 8.0 for those objects. 

The temperatures measured from the SDSS photome- 
try were found to be systematically lower than those from 
the fitting of line profiles (e.g. Table[9] objects with SDSS 
spectroscopy, but fitted with the photometric method). Our 



sample of SDSS white dwarfs overlaps with the Palomar 
Green (PG) sample, and Fig.[7]shows a comparison between 
ou r spectroscopi c and photometric temperatures with those 
of iLiebert et~atl (|2005l ). which were determined from inde- 
pendent data, spectral models, and fitt ing routines. We find 
good agreement between the results of ILiebert et al.l (|2005l ) 
and our spectroscopic method. However, our photometric 
temperatures are systematically too low, a trend that is 
strongly correlated with either white dwarf temperature or 
distance. At 200pc (500pc), the photometric temperatures 
are on average 5% (10%) below our and Liebert's spectro- 
scopic values. This could suggest that interstellar redden- 
ing is, at least in part, the culprit for reduced tempera- 
tures. While reddening would not significantly effect the 
shape of the line profiles, it cou ld noticeably chang e the 
slope of the continuum (see also iHolberg et all l200Sl ). No 
clear correlation is, however, s een when comparing t he mis- 
match in temperatures to the ISchlegel et al.l l|l998h values 
of E(B — V) at the positions of the white dwarfs. The 
ISchlegel et al.l maps probe interstellar reddening through the 
entire Galaxy, whereas the white dwarfs in our sample lie 
at a few hundred parsecs at most. Typical (total Schlegel) 
reddening along the lines-of-sight towards our white dwarfs 
is E(B - V) ~ 0.05. De-reddening the SDSS photome- 
try with that total E(B — V), and re-fitting the photo- 
metric white dwarf sample indeed leads, as expected, to 
a large over-correction of the white dwarf temperatures. 
Analysing the sample of spectroscopic DAs, we estimate 
that the typical reddening in front of the white dwarfs is 
E(B — V) ~ 0.01 — 0.02. However, we can not systemat- 
ically correct for the effect of reddening for the sample of 
photometric-only DA candidates. We note that for hot white 
dwarfs, T e ff > 45000K, non-LTE effects become important, 
which may also lead to some systematic differences in the fit 
parameter for the hottest stars in our sample. Therefore, the 
temperatures calculated from the photometry alone have an 
additional systematic uncertainty, on top of the statistical 
uncertainty from the fit, and the true temperatures are likely 
to be a few thousand Kelvin higher. In the context of our 
search for infrared flux excess (Sect. [5721 below), changing 
the white dwarf temperature by a few thousand Kelvin does 
not have a significant impact on the level of excess detected 
(see Sect. [6] for examples). 



5.2 IR Excess Detection 

In both the spectroscopic and photometric fitting methods, 
the best fit model extends into the IR and all objects with 
UKIDSS data were examined for an excess by comparing 
the IR white dwarf model flux with the observed YJHK S 
magnitudes. Objects with a 3<r excess in any band over the 
white dwarf model were defined as a robust excess candidate 
("DAire" and "DA:ire" for the spectroscopically confirmed 
DAs and the photometric-only DA candidates, see Tables [7J 
[8] and [9}. Further to this, objects that appeared to have a 
best fit model which over-estimated the flux in Y, J and H, 
but showed only a ~ 2cr excess in K were also flagged as ten- 
tative excess candidates ("DAire:" and "DA:ire:", as above, 
Tables[7l|8]and|9}, see Fig.[5]for an example. Further IR data 
is definitely needed to confirm these marginal IR excess can- 
didates. Similarly, spectroscopically confirmed DAs (photo- 
metric DA candidates) with close to 3cr excess that by eye 
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Figure 6. SDSS J0135+1445; a cool white dwarf with a probable low mass companion. The fit is a T cff = 7467±18K, log(g) = 7.34±0.04 
spectra at a distance of 69 ± 2pc. The calculated mass is M w( j = 0.29 ± O.OIMq. The photometry is best fit with a white dwarf of 
T e ff = 8000±jqK, log(g) = 8.0. Based upon the fit to the photometry, the addition of an L8-type companion to the white dwarf 
magnitudes is the best fit to the IR excess (Table[9]| and is plotted on the right hand figure as open grey circles. The model Y-band 
magnitude does not match the UKIDSS Y measurement, but we found that the Y-band data did relatively often disagree with SDSS z 
and UKIDSS J. 




PG T eff [1000K] Spectroscopic d [kpc] 

Fig ure 7. A demonstra tion of the systematically lower measured photometric effective temperatures. The left hand panel shows 
the iLiebert et al. PG survey wh ite dwarfs with corr esponding stars in our sample (25 objects). This figure compares 

the effective temperature calculated by ILiebert et al.l l l2005t) with that from our spectroscopic (grey) and photometric (blue) 
fitting techniques. When using the photometric method, we significantly underestimate the temperature of hot white dwarfs. The 
hottest white dwarf at 70, 000K has a large error on the photometric fit because the continuum slope of the white dwarf at these 
temperatures is effectively Rayleigh-Jeans and therefore contains no temperature information. In the right hand panel we show 
the ratio of the effective temperatures calculated using the photometric and spectroscopic fitting methods for each of the objects 
in the spectroscopic DA white dwarf sample, as a function of distance. Again, we underestimate the effective temperature of white 
dwarfs at large distances. 



require further data to confirm the excess were also marked 
as "DAire:" ("DA:ire:"). For the photometric-only objects 
the uncertainty on the model parameters is generally larger 
compared to the spectroscopically confirmed DAs. This was 
accounted for by not flagging objects with a marginal IR 
excess and a large uncertainty on effective temperature as 
excess candidates. 



Spurious excesses were often caused by spatially close 
background or foreground objects to the white dwarf and 
bad SDSS or UKIDSS images. We visually inspected all 
flagged sources and discounted resolved and partially re- 
solved systems because the physical association of the two 
objects could not be demonstrated based on the available 
data. 
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A total of 42 white dwarfs were found to have an ex- 
cess from the spectroscopic fitting method, and 105 infrared 
excess candidates were found from the photometric fitting 
method (Tabled . The excesses have a variety of spectral 
shapes, generally consistent with various spectral type com- 
panions. Table[7] and [8] list all the spectroscopic and photo- 
metric IR excess candidates, respectively. 

5.3 IR Excess Modelling 

The YJHK„ magnitudes of objects that were found to ex- 
hibit an infrared excess were fitted with a composite model 
consisting of the best-fit white dwarf plus a set of low-mass 
companions with spectral types MO through to L8. W e used 
the 2MASS JHK S magnitudes of lHoard et aD (j2007h . which 
we converted into the UKIDSS fil t er system ado pting the 
equations given in iDve et all (|2006h ■ iHoard et all s absolute 
magnitudes of the low-mass companions were scaled to ap- 
parent magnitudes using the distance modulus calculated 
from the white dwarf fit. For photometrically fitted objects, 
the white dwarf distance modulus was calculated as an av- 
erage of the difference between the best-fit absolute mag- 
nitudes and the SDSS apparent magnitudes in each of the 
u, g, r and i bands. Finally, the model magnitudes of the 
composite system were computed from the combined white 
dwarf and companion star fluxes. The best composite fit 
to the JHK IR photometry was calculated using a least 
X 2 search and subsequently confirmed by visual inspection. 
Fig|8] shows reduced \ 2 as a function of companion type for 
three white dwarfs with well constrained companions using 
the photometric method as an example. It was seen that a 
white dwarf with a good fit to the companion has a reduced 
X 2 10, an example of a the corresponding SED is shown 
in Fig. [6] Composite fits with \ 2 ~ 10 — 100 were flagged 
as "bad fits" and flagged as DAire: or DA:ire in Table[7] 
and |SJ respectively, as the nature of the IR excess remains 
somewhat unclear (see Fig.|9]and Sect . 15. 4T) . Finally, QSOs 
stand out because of their very high X 2 , ~ 100 - 1000 (e.g. 
Fig. I10p . A summary of companion types for each method is 
given in Table[4] 

5.4 Comparison of the spectroscopic and 
photometric methods 

The results from using either the spectroscopic or purely 
photometric methods to find IR excess can be directly com- 
pared because all the objects with SDSS spectra also have 
SDSS photometry. We discuss in the following sub-sections 
how well the spectroscopic and photometric method com- 
pares for genuine DA white dwarfs (Sect. 15.4.1)1 . what can 
be learned from the spectroscopic contaminants that show 
an IR excess when fitted with the photometric method 
fSect. l5~4~2l 15,4.3"]) . and how additional information such as 
IR colours and proper motion could be used to suppress 
contaminants in the photometric-only sample fSect. [574~2j) . 

5.4.1 DA White Dwarfs 

Of the 42 spectroscopic DA white dwarf IR excess candi- 
dates, 30 (~ 71%) are recovered by the photometric fitting 
method as well (column "DA w SIRE" in Tableg] 9 M-type, 
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Figure 8. Reduced \ 2 as a function of companion type for three 
white dwarfs with well constrained companions using the photo- 
metric method. \ 2 ls calculated from comparing the UKIDSS 
JHK magn it udes with those of the low-mass companions of 
IHoard et all 11200711 . SDSS J0135+1145, SDSS J0842+0004 and 
SDSS J0925-0140 are plotted as blue, grey and black respectively. 
We find that a good fit has a reduced \ 2 < 10, a bad fit (DAire: 
or DA:ire: in Tables[7] and [Sj has \ 2 ~ 10 - 100 and a QSO has 
X 2 ~ 100 - 1000. 
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Figure 9. The SED of SDSS J1619+2533, an example of a spec- 
troscopically confirmed DA white dwarf where the photometric 
method substantially underestimates the white dwarf temper- 
ature. The best fit to the spectrum is a T e ff = 25595K and 
log(g) = 7.21 DA white dwarf model (M wd = 0.33 ± 0.04), 
where as the photometric fitting method finds the best solution 
at T e fj = 18000K and log(g) = 9.5, the latter is plotted here in 
black. Within the photometric method, the lower temperature 
(and higher gravity) leads to the distance being substantially 
underestimated, and therefore the flux of the companion being 
overestimated. Fitting the SED with the two-component model 
described in Sect. 15.31 compensates for the low distance by choos- 
ing a companion with a larger absolute magnitude, i.e. a later 
spectral type. For SDSS J1619+2533, the spectroscopic method 
results in a M4 type companion, the photometric method in an 
L6 companion. Because the optical-IR SED can not be well fitted 
with any companion type at the photometric distance, this object 
is flagged as "bad fit" in Table [9] 
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Table 4. Objects with an IR excess split by estimated low-mass companion type. The first column shows the IR excess candidates among 
the spectroscopically confirmed DAs (Table[7Jl and the final column contains the information on the photometric-only sample (Table[HJl, 
where there is no SDSS spectrum available for classification. The section in the middle shows all objects with SDSS spectroscopy for 
which fitting their ugri photometry resulted in an IR flux excess (Table[9]l, split by their spectroscopic classification. This set of objects is 
useful to gauge the contamination of photometric-only IR excess DA candidates by other types of objects. It also allows the like-for-likc 
comparison of the infrared excess detection/modelling for genuine DA white dwarfs. To aid this comparison of the two methods, the 
column "DA (w SIRE)" lists the spectroscopically confirmed DA white dwarfs that exhibit an IR flux excess when fitted with both the 
spectroscopic and photometric method, and the column "DA (wo SIRE)" lists 6 spectroscopically confirmed DA white dwarfs that are 
found to have an IR excess only when analysed with the photometric method. A "-" mark denotes that the corresponding entry is not 
possible. "No Phot IR Excess" means that no photometric excess was found, even though a spectroscopic excess is seen. 



Spectroscopic Spectroscopic using Photometric 

Photometric Method 



Primary Type => 
Companion Type ij. 


DA 


DA 

w SIRE wo SIRE 


NLHS 


WDMS 


DAH 


QSO 


Total 




M-type 


16 


9 


1 


10 


2 








22 


15 


L-type 


19 


11 





12 


1 








21 


19 


^ L7 or disc 


7 


6 


2 














8 


4 


QSO 



















323 


323 


40 


Bad fit 





1 


3 


9 


1 


3 





20 


27 


No Phot IR Excess 




12 












12 




S 


42 


42 


6 


31 


4 


3 


323 


409 


105 



11 L-type, 6 ^ L 7- type companions/debris discs and 4 "bad 
fits", which are marked as "DAire:" in Table[8]l, and 12 do 
not show an excess when using the photometric method. 
These 12 objects are close to the ~ 3a limit in the spec- 
troscopic method and were not recovered. They were not 
recovered by the photometric method because of the larger 
uncertainty on effective temperature leading to the excesses 
being within the combined model and flux errors. 6 objects 
with DA white dwarf SDSS spectra are found to have an IR 
excess in the photometric method, but do not exhibit an IR 
excess when using the spectroscopic method (column "DA 
wo SIRE" in TableU) . Three of these have excesses which are 
accentuated by a slightly hotter photometric fit compared 
to the spectroscopic one. They fall just outside the crite- 
ria (Section l5.2|) for having an excess in the spectroscopic 
method. For the other three objects the SDSS spectra are 
too poor to obtain reliable T c ff and log(g) measurements 
from fitting their Balmer lines. Therefore an excess was not 
recognised in the spectroscopic method. However, the photo- 
metric data were good enough and thus these three objects 
are classified as photometric IR excess candidates. We may 
therefore expect ~ 9 — 17% (3/36-6/36) spurious IR excess 
candidates from the photometric method. All objects that 
display an infrared excess either in the spectroscopic, or the 
photometric method, and have an SDSS spectrum are listed 
in Tabled 

Table[4] also lists the IR excess candidates split by mod- 
elled companion type. Among the 30 objects which are de- 
fined as IR excess candidates in both methods, the distri- 
bution of companion type has a similar form. ~ 80% of the 
companions have M or early L spectral types (split evenly 
between the two classes), and the remaining ~ 20% have 
companions ^ L7, i.e. brown dwarf or debris disc candi- 
dates. 

As briefly outlined in Sect . 15.11 the temperatures result- 
ing from the photometric fits are systematically too low, 
and this will introduce a bias in the spectral type of the 



companion. The flatter SED of a cooler white dwarf will re- 
duce/distort the IR flux excess relative to the white dwarf, 
and therefore a later type companion will provide sufficient 
flux to account for the excess. An additional effect is that 
a cooler white dwarf will suggest a lower distance when fit- 
ting the ugri magnitudes. Underestimating the distance will 
lead to an overestimate of the absolute flux of the compan- 
ion. To compensate for this, the fit to the companion will 
resort to a companion with a larger absolute magnitude, i.e. 
lead to a companion spectral type that is too late. A mod- 
erately extreme example of these effects is shown in Fig.[9] 
(SDSS J1619+2533, see Table©. Assuming that the spec- 
troscopic fit parameters are correct, the white dwarf tem- 
perature is underestimated by ~ 8000K, which leads to an 
L6 companion in the photometric method, as opposed to 
an M6 companion resulting from the spectroscopic analysis. 
This object is marked in the notes column as having a bad 
fit to the companion in the photometric method, where a 
"bad fit" is defined as an object having an excess which is 
inconsistent with any companion type (at the photometric 
white dwarf distance). 

5.4.2 Quasar Elimination 

As can be seen in Table[4] a significant fraction (323) of 
the IR excess candidates from the photometric method with 
spectra are quasars. This is almost the entire population of 
spectroscopically confirmed quasars in the DA white dwarf 
sample with UKIDSS magnitudes. The remaining 5 quasars 
do not have sufficient IR data to show their quasar nature, 
but equally are not flagged as IR excess can didates. This 
is cau sed by the flat SED exhibited by QSOs (I Covey et al.l 
120071 ). When modelled with a DA white dwarf, they gener- 
ally have an effective temperature of ~ 8000K and have an 
IR excess that is much higher than physically possible for an 
M-dwarf or brown dwarf companion at the distance of the 
(photometric) DA fit (e.g. Fig. llOf) . Because of the (appar- 
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Figure 10. The SED of SDSS J0046-0044, an example of a 
quasar, that was selected as a candidate DA white dwarf (based 
on its colours, Table[T} with possible IR excess ("DA:ire:") by 
the photometric method. The best-fit to the ugri photometry is 
found for T c g = 9000K and log(g) = 7.0 at a distance of 37pc 
and shown as black line. The excess over the model keeps rising 
steeply into the mid-IR and can not be modelled by any compan- 
ion type at the distance of the (photometric) DA fit, identifying 
this object as a quasar contaminant. 

ent) low effective temperatures, the quasars also generally 
have low (apparent) distances. These properties, however, 
overlap with those of genuine cool nearby white dwarfs and 
cannot be used to directly distinguish between white dwarfs 
and quasar. Another common sign of QSOs are large jumps 
in brightness between adjacent magnitudes (caused by emis- 
sion lines) that are not seen in any of the genuine DA white 
dwarfs with infrared excess. Given that we correctly clas- 
sified, on the base of their optical-IR SED, 99% of (spec- 
troscopic) quasars contaminants that were found as DA:ire 
by the photometric method, we are confident that we can 
identify the vast majority of quasars among the photometric- 
only sample. 

Among the photometric-only sample of IR excess can- 
didates, we find 38 objects whose SEDs are very similar to 
those of our 323-strong spectroscopicially confirmed quasar 
sample and which we therefore believe to be quasars as 
well. Optical spectroscopy is needed to confirm their nature. 
These have been removed from Table[8]and can be found in 
the online QSO tables available via CDS. 

5.4-3 Contamination by NLHS and Non-DA White 
Dwarfs 

Contaminant NLHS and non-DA white dwarfs are more 
difficult to identify and remove from the photometric-only 
DA candidate sample, as their overall SEDs are all rather 
similar. Closer inspection of the SDSS spectroscopy of the 
NLHS and cross-checking them in Simbad suggests that a 
large fraction of them are subdwarf B (sdB) stars. The bulk 
of subd warfs are believe d to have for med in binary inter- 
actions (|Han et all 120031 ; iHeberl |2009l ) and therefore it is 
expected that a large majority will still have companions. 
Such companions would cause an IR flux excess over the 
Rayleigh- Jeans tail of the subdwarf, and it is hence likely 



that our photometric-only sample of candidate DA white 
dwarfs with infrared excess contains a significant contami- 
nation from sdB plus low-mass companion binaries. 

Given that we fit all photometric objects with DA model 
spectra, we may expect some imperfections in the fits to the 
photometric NLHS objects. Nevertheless, the detection of a 
near-IR flux excess over the Rayleigh- Jeans tail of the model 
is likely to be correct for many of the NLHS objects among 
the photometric sample. However, a DA fit to the photom- 
etry of a physically much larger NLHS object will dramati- 
cally underestimate its radius, and hence its distance. Con- 
sequently, fitting the companion with the composite model 
(Sect. [573)) will result in a spectral type of the companion 
that is much too late. 

From the spectroscopic sample, we find that the fre- 
quency of NLHS inside the DA colour selection (Table[T] 
Fig-IUl is 14.6% (TablefSJ, however the level of contamina- 
tion is strongly magnitude-dependent. Figure[TT](left) shows 
the distribution of NLHS and DA white dwarfs as a func- 
tion of (?-band magnitude. The ratio of these two, and thus 
the expected level of contamination of the photometric-only 
sample, is shown in the middle panel (this assumes that 
the majority of quasar have been removed because of the 
characteristic shape of their SED). The contamination of 
the photometric-only DA candidate sample by NLHS drops 
significantly towards fainter g magnitudes. Subdwarfs are 
~ 100 times brighter than white dwarfs and therefore ap- 
parent magnitudes of g ~ 18 — 19 sample distances of many 
kpc. This is several times the exponential scale height of 
the Galactic thick-disc population, and hence the number 
of sdBs at such large distances is relatively small. Figure[TT] 
(right panel) shows the photometric-only IR excess candi- 
dates as a function of g magnitude, where we can assume 
that most objects with g < 16 are likely to be NLHS. 

An additional clue on the NLHS vs DA classification 
of the photometric-only objects comes from their location 
in the (u — g,g — r) colour-colour diagram (Fig. 1121 bottom 
right panel), where the majority of the spectroscopically con- 
firmed NLHS objects are concentrated at the blue end of the 
"DA" banana. We can therefore assume that NLHS are the 
primary contaminants of the photometric sample in this re- 
gion as well (Fig. 1121 top right panel). 

Assuming that the distribution of NLHS contaminants 
is similar between the spectroscopic and photometric-only 
sample, we expect ~ 259 of the 1771 photometric objects to 
be NLHS. When fitting the photometric sample, 31 (12%) 
spectroscopically confirmed NLHS were found to have an IR 
excess (red dots in the bottom right panel of Fig. ll2[) . Again, 
assuming that the contamination among the spectroscopic 
and photometric samples is similar, we would expect that 
~ 37 (12% of 259) of our photometric DA candidates with in- 
frared excess (DA:ire and DA:ire:) are in reality NLHS, pri- 
marily sdB stars with main sequence companions. These are 
still interesting in their own right (e.g. Sect. l6.3.2.1|l . but not 
the primary focus of the IR excess search (a detailed inves- 
tigation of subdwarfs with main-sequence star companions 
will be published elsewhere, Girven et al. (2011, in prep). 

In contrast to the above discussion on NLHS, fitting 
non-DA white dwarfs with DA white dwarf models is likely 
to provide a reasonably good estimate of all system param- 
eters, including the companion type - the downside being 
that on the base of photometry only, it is nearly impossible 
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to differentiate between DA and no n-DA white dwarf s . How - 
ever, based on the statistics of the lEisenstein et al.l l|200rj ). 
we only expect a small level of contamination by non-DA 
white dwarfs (Table[2j). 



5.4-4 Independent checks: proper motions and infrared 
colours 

Our method of identifying infra- red excess candidates fol- 
lows [lYembla^&!Bergeron[ ([2007]) , i.e. fitting model spectra 
to SDSS spectroscopy and photometry. In this section, we 
carry out an independent investigation of our sample using 
colour-colour diagram s, s uch as previo u sly ex plored by e.g. 
IWachter et all i|2003l ) and lHoard et al l (2007), and test our 
classification of the photometric-only systems by making use 
of proper motions. 

Figure[l3] shows the distribution of the spectroscopic 
SDSS/UKIDSS sample in the (z-H,H- K) colour-colour 
space (left panel). Using model white dwarf colours, it can 
be seen that the DA cooling sequence runs from top left 
to bottom right through the white dwarf group. The NLHS 
and DA white dwarfs are clearly separated from the quasars, 
which can be understood as stars are on the Rayleigh- Jeans 
tail i n the infrared, wh ereas quasar follow a flatter power 
law j|Covev et al.l 120071 ). Therefore, as an additional test 
of how reliable our classification of the photometric-only 
quasar candidates works, we inspected the infrared colours 
of the photometric-only sample. Choosing an empirical cut 
of H — K > 0.627 selects 89.6% of the quasar contaminants 
in the spectroscopic sample. Adopting the same H — K cut 
for the photometric-only sample flags 38 objects as quasar 
candidates. This includes 31 of the 38 photometric-only ob- 
jects identified as quasar candidates on the base of their 
SED (Sect. l5A2l (see online QSO table), of which only 37 
have both H and K measurements). In addition, two of 
the H — K selected QSO candidates correspond to the two 
"weak" photometric-only quasar candidates listed in Ta- 
ble! 

A third concentration of objects can be seen in Fig. 1131 
in the region between the DA white dwarfs and the quasars, 
a significant fraction of which are DA infrared excess can- 
didates. The presence of an infrared excess over the stellar 
flux distribution results in a flatter spectral slope, and hence 
moves these objects closer to the quasar locus. On closer in- 
spection, some white dwarfs in this region are found to be 
blended sources or suffer from background contamination 
from a nearby galaxy and therefore were not included as 
infrared excess candidates. 

In summary, adopting H — K > 0.627 would effi- 
ciently remove the bulk of contaminating quasars from the 
photometric-only sample, however, such a cut would also 
remove a handful of genuine white dwarfs with the largest 
infrared excess emission (such as e.g. SDSS J1228+1040, see 
Sect. ltTXTj) . 

Given that white dwarfs are nearby low-luminosity ob- 
jects, they are expected to exhibit larger proper motions 
than the more luminous NLHS, and quasars are not ex- 
pected to show any significant proper motion at all. Thus, 
proper motions can be used to distinguish between white 
dwarfs and quasars in the region with H — K > 0.627. We 
have retrieved proper motions from DR7 for all objects in 
the SDSS/UKIDSS sample. Figure[TJ (left panel) shows a 



cumulative proper motion distribution for the white dwarfs, 
quasars and NLHSs. Based upon this, we chose a cut in 
proper motion of p.m. ^ lOmas/yr to define low proper 
motion objects such as quasars. In the spectroscopic sample, 
this cut selects 7%, 74% and 97% of the DA white dwarfs, 
NLHS and quasar respectively, efficiently eliminating the 
majority of the quasar without removing too many white 
dwarfs. 

FigureQT] (middle panel) plots the magnitude of the 
proper motion as a function of H — K. The statistical sig- 
nificance of the proper motion is encoded in the size of the 
points, where larger points denote more significant proper 
motions. As expected, the spectroscopically objects classi- 
fied as NLHS stars and quasars show very small proper 
motions, which are in most cases consistent with zero. The 
right hand side panel of Fig.[l4] shows the location of 35 
photometric-only objects classified as quasar because of the 
characteristic shape of their SEDs (see online QSO ta- 
ble) that have both H — K colours and proper motions. 
The vast majority of these objects are contained within 
H — K > 0.627 and p.m. ^ lOmas/yr, corroborating our 
SED-based classification. 

A small number of objects with H — K > 0.627, i.e. 
within the "quasar" region, display large and statistically 
significant proper motions (Fig. 1141 middle panel). These are 
listed in Table[5] Among those objects are three spectroscop- 
ically confirmed DA white dwarfs. At closer inspection, the 
UKIDSS magnitudes of SDSS J1244+0402 may be contam- 
inated by a nearby background object, which would lead 
to a spurious H — K colour. The other two spectroscopic 
DA white dwarfs, SDSS J0753+2447 and SDSS J1557+0916 
have very red H — K colours and high proper motions. 
They are therefore excellent IR excess candidates. Apply- 
ing the same procedure to the photometric-only DA can- 
didate sample, SDSS J0959— 0200 is the strongest infrared 
excess candidate among the three photometric-only objects. 
SDSS J1440+1223 and SDSS J1509+0539 also appear to suf- 
fer from background contamination. 

Another four spectroscopic objects with quasar-like IR 
colours and high proper motions are classified as three mag- 
netic white dwarfs (DAH) and one DZ white dwarf. We 
would expect that the IR spectra of all of these objects 
should be close to a Rayleigh- Jeans distribution, suggest- 
ing that the inferred IR excess is probably real. In fact, one 
of the DAH, SDSS J 1212+0136 is a well- studied DAH plus 
brown dwarf binary dSchmidt et al.ll2005Tl . which exhibits a 
genuine infrared excess l|Debes et al.ll2006l l (see Sect. RTO)) . 
As discussed in the introduction, the metals seen in DAZ 
white dw arf atmospheres are from recent or ongoing accre- 
tion (e.g. iDupuis et al.lll993l ; iKoester fc Wilkenll2006h. The 
same is tr ue for cool DZ white dwarfs (e.g. lAannestad et"al] 
Il993l ). and lFarihi et"aH l|20ld ) discuss the potential connec- 
tion between the large number of DZ white dwarfs, and the 
DAZ white dwarfs with dusty debris discs thought to orig- 
inate from the tidal disruption of rocky asteroids. The DZ 
found here, SDSS J1342+0522, exhibits a very red H - K 
colour. The excess in K over the model spectrum is slightly 
under 3<r and therefore is, in our classification scheme, only 
a marginal candidate for having an IR excess. This DZ white 
dwarf and the two new DAH white dwarfs IR excess candi- 
dates warrant further investigation. 

Finally, eight quasars have statistically significant 
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Figure 11. Left panel: The distribution of DA white dwarfs (blue, left hand scale) and NLHS (black, right hand scale) as a function 
of (/-band magnitude. Middle panel: the fraction of NLHS as a function of <?-band magnitude. The probability of contamination 
by NLHS clearly drops with increasing g-band magnitude. Right panel: distribution of the photometric-only IR excess candidates, 
the brightest of these group are most likely NLHS rather than DA white dwarfs. 



proper motions (~ 3 — 5a; see TabreJSJ, which highlights 
the fact that the SDSS vs USNO-B proper motions have 
to be considered with a pinch of salt: among a total of 328 
quasars with SDSS spectra and UKIDSS data, we would ex- 
pect only one to have a 3a significant proper motion, and 
none at 4a. 

We conclude that dissecting the white dwarf sample se- 
lected with our constraint set (Table[l]) using colour-colour 
diagrams and proper motions leads to mutually consis- 
tent results when compared to our primary methodology 
(Sect. [5711 - 15. 3|) . but the spectroscopic modelling provides an 
additional wealth of information. 



5.5 Overall Numbers 

A summary of the numbers of objects at each stage of the 
analysis is given in Table[6] The total numbers of white 
dwarfs with near-IR excess are broken down according to 
the spectral type of their companions in Table|4] 

We find that 3.3% (42 of 1275) of the SDSS spectro- 
scopically confirmed DA white dwarfs with at least one of 
the YJHK UKIDSS magnitudes have an IR excess and are 
therefore candidates for having a companion or a debris disc. 
However, this does not take account of the fact that the sam- 
ple of white dwarfs is incomplete even within UKIDSS DR8 
because for many only subsets of the IR magnitudes are 
available. Thus we are limited by UKIDSS coverage and the 
real number is higher. 2.0% of the spectroscopic DA white 
dwarfs are candidates for having a companion of type L0 or 
later, i.e. have brown dwarf companions. Similarly, 0.5% are 
promising disc candidates, having an excess compatible with 
a companion type of L7 or later. If we only discuss the ob- 
jects with a detection in the K-band (required for detecting 
a disk), 1.2% are disk candidates, and only including objects 
where we are confident of the IR excess (not "DAire:" in Ta- 
ble0 see Sectionl531 and \5~3\i , a lower limit of 0.8% of DA 
white dwarfs have a brown dwarf companion. 

For the photometric-only sample, where we have fit- 
ted the SDSS and UKIDSS photometry, 5.9% are IR ex- 
cess candidates. However, this number will be affected by 
the efficiency of white dwarf selection (62.3%) and the effi- 
ciency of removing contaminants. Assuming that we remove 



Table 5. Objects with quasar-like infrared colours (H — K ^ 
0.627), but large lOmas/yr) and statistically significant (> 3ct) 
proper motions. The significance of the proper motions is listed 
as CTp.m.- Classifications are given for the objects that have SDSS 
spectra. The comment "BG object" refers to a second resolved 
source being seen in the UKIDSS images. These are most likely 
background galaxies which would significantly affect the H — K 
colour. For full coordinates, refer to the online spectroscopic and 
photometric-only tables available via CDS. 



Name 


H - K 


p.m. (mas/yr) 




Class 


Comment 


0043+0005 


0.79 ± 0.06 


13.68 ± 3.06 


4.5 


QSO 




0753+2447 


0.64 ± 0.24 


32.16 + 3.20 


10.0 


DA 




0858+0938 


0.88 ± 0.03 


13.13 + 3.01 


4.4 


QSO 




0959-0200 


0.65 ± 0.15 


30.25 ± 3.10 


9.8 






1031+0341 


1.05 + 0.08 


12.25 ± 3.92 


3.1 


QSO 




1142+1347 


1.02 ± 0.05 


10.10 + 2.81 


3.6 


QSO 




1212+0136 


1.00 + 0.06 


69.19 + 3.21 


21.6 


DAH 


(1) 


1244+0402 


0.88 + 0.17 


11.80 + 3.09 


3.8 


DA 


BG object 


1250+1549 


1.49 ± 0.02 


93.14 + 3.12 


29.9 


DAH 




1342+0522 


0.64 + 0.27 


51.53 + 3.24 


15.9 


DZ 




1427-0054 


0.68 ± 0.04 


15.99 + 3.25 


4.9 


QSO 




1440+1223 


1.43 ± 0.08 


88.83 + 2.65 


33.5 




BG object 


1443+0910 


0.96 ± 0.02 


12.46 + 2.98 


4.2 


QSO 




1509+0539 


0.98 + 0.25 


48.29 + 3.07 


15.7 




BG object 


1514+0744 


0.76 + 0.07 


102.64 + 5.06 


20.3 


DAH 




1553+0718 


0.80 + 0.01 


14.64 + 2.75 


5.3 


QSO 




1557+0916 


0.70 + 0.21 


25.48 + 2.85 


8.9 


DA 




1557+2646 


0.84 ± 0.04 


10.07 + 3.21 


3.1 


QSO 





(1) A magnetic cataclysmic variable in a low s t ate llSchmidt et all 
20051: iDebes et al.ll2006l; iBurleigh et al.ll2006l ; iFarihi et al.ll2008l: 
Howell et alj|2008l : iLinnell et alj|20ld) . 



all the obvious photometric-only quasar contaminants (38, 
see online QSO table) from the IR excess objects, the re- 
maining 67 photometric-only IR excess candidates will be 
either DA white dwarfs or NLHS. As discussed previously, 
we estimate that ~ 37 of these are NLHS (Section l5.4.3p . 
This number is very similar to the amount of "bad fits" 
found in from the photometric method and so we believe 
the majority of "bad fits" come about from NLHS contami- 
nants in the photometric-only sample. Therefore, we expect 
~ 30 (2.7% of ~ 1103) genuine DA white dwarfs with IR 
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Figure 12. Top left panel: the distribution of the 42 spectroscopically confirmed DA white dwarfs with infrared excess (Table[7jl 
in the (u — g, g — r) colour-colour space. The black dashed line outlines the (u — g,g — r) DA colour selection (Table[T] ig.[2j| . Top 
right panel: the distribution of the 74 IR excess objects found by photometric method (excluding quasars and quasar candidates), 
with DA white dwarfs, non-DA white dwarfs and NLHS shown in blue, yellow, and black respectively. SDSS J1228+1040 and 
SDSS J1043+0855, two DA white dwarfs known to have debris discs, are indicated by pink and red bulls-eye symbol, respectively 
(top left and right panel, see Sect. l6.1~Tl 16.1.21 . SDSS J1212+0136, a known magnetic white dwarf with a brown dwarf companion 
is shown as purple bulls-eye symbol (right panel only, see Sect, l6.1~3)l . Bottom left panel: the distribution of the 67 photometric- 
only infrared excess candidates. Bottom right panel: all photometric-only objects satisfying our DA constraint set from Table[T] 
(1771 objects, gray dots), spectroscopically classified NLHS (209 objects, black), and spectroscopically classified NLHS with an 
infrared excess (31 objects, red). Based on the analysis of the spectroscopic sample (Sect. [574.31 1. we estimate that ~ 12% of the 
photometric-only objects are NLHS, with a strong concentration towards the (blue) top-end of the DA "banana". 



excess among the photometric-only DA white dwarf candi- 
dates. These infrared excess candidates show a similar dis- 
tribution in companion type when compared to the spectro- 
scopic sample. Considering only the objects we are confident 
of the excess (not "DA:ire:" in Table[5J see Section [5T2l and 
15. 3|) . 1.8% (19) of the photometric-only DA white dwarfs 
candidates are likely to have a brown dwarf companion. 



6 NOTES ON INDIVIDUAL OBJECTS 

Some objects of particular interest identified in the SDSS / 
UKIDSS cross-correlation are discussed below. We separate 
those into objects already known to host a disc or companion 
(Section l6.1[) . notable objects from the spectroscopic sample 
(Section l6.2l) . and the same from the photometric-only sam- 
ple (Section l6.3[) . 



6.1 Prototypical White Dwarfs 

Our spectroscopic sample contains two white dwarfs that 
were known to host debris discs and one magnetic white 
dwarf with a substellar companion, and they serve hence as 
a benchmark for our selection procedures. 



6.1.1 SDSSJ1228+1040 

SDSS J1228+1040 is one of two DA white dwarfs in our 
SDSS/UKIDSS sample known to have a debris disc. The 
disc was initially identified because of the highly unusual 
emissi on lines of the Ca II 8200 A triplet (Ga nsicke et al.l 
l2006bl ). whose double-peaked shaped can only be explained 
by metal-rich gas orbiting the white dwa rf within its tidal 
disruption radius ([Gansicke et alJ l2006bl ). Near- and mid- 
IR observations revealed a substantial infrared excess over 
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Figure 13. Location of the SDSS/UKIDSS sample in (z — H, H — K) colour space. The left panel shows the objects with SDSS 
spectra, DA white dwarfs, NLHS, other white dwarfs and quasars in blue, black, cyan and green, respectively. The 23 robust 
"DAire" infrared excess objects from Table[7]are shown as red open circles. The positions of, SDSS J1228+1040, SDSS J1043+0855 
and SDSS J1212+0136, two white dwarfs in the sample known to have a debris disc and a DAH white dwarf plus brown dwarf 
binary (see Sect I6.1TT1 16.1.21 16.1.3j l. are indicated by pink, red and purple bulls-eye symbols, respectively. A track of model DA 
white dwarf colours for 6000K T e fj ^ 100000K and log(</) = 8.0 is plotted as a solid black line. High temperature white dwarfs 
are found at the top left end of the line. The scatter about this line cannot be explained by a spread in log(g) alone, but by 
uncertainties on the magnitudes (primarily those from UKIDSS) and background contamination. The photometric-only objects 
are shown in the right hand panel, with the 29 robust "DA:ire" infrared excess candidates (Table[8j) again highlighted in red. 




PM [mas/yr] H-K H-K 



Figure 14. Left panel: The cumulative histogram of spectroscopically confirmed DA white dwarfs (blue), NLHS (black) and 
quasar (green) as a function of proper motion (PM). Black dashed lines show a proper motion of lOmas/yr and corresponding 
population levels. 7% of DA white dwarfs are contained within the bin of PM ^ lOmas/yr. Similarly, 74% of NLHS and and 
97% of quasar. Middle panel: distribution of proper motions as a function of H — K colour, using the same colour-coding, and 
showing in addition non-DA white dwarfs in cyan. The significance of the proper motion is encoded in the size of the points, where 
larger points denote more significant proper motions, clipped at a maximum of 30ct. Black dashed lines show PM = lOmas/yr 
and H ~ K = 0.627. The position of the three benchmark objects SDSS J1228+1040, SDSS J1043+0855 and SDSS J1212+0136 
(see Sect 16, lTTl 16.1,21 16. 1,31) are indicated by pink, red and purple bulls-eye symbols respectively. Right panel: same as the middle 
panel, but showing the (35) quasar candidates identified on the basis of their SED shape (see online QSO table) which have both 
H and K band data in UKIDSS and proper motions in SDSS DR7. 



the white dwarf, unambiguously identifying a dusty com- 
ponent of the debris d isc, in addition to the gaseous one 
ijBrinkworth et al"1l2009l ). Our fits to the SDSS spectroscopy 
and photometry are shown in Fig. 1151 The two fits differ in 
T e ff by 2000 K. However, this has very little effect on the 
extrapolated infrared magnitudes of the white dwarf, and 
the object shows a 3cr excess in H and 12a in K, indepen- 
dent of the method used for fitting. At the temperature of 
SDSS J1228+1040, \og{g) is not well constrained from fitting 
the ugri photometry alone, however, this primarily effects 



the widths of the hydrogen lines, and has a negligible ef- 
fect on the spectral slope of the white dwarf model. This 
demonstrates that the detection of a genuine infrared excess 
is robust and independent of whether an optical spectrum 
is available. 

Taking the infrared measurements at face value and 
ignoring our knowledge about this star, we have modelled 
the SDSS/UKIDSS spectral energy distribution as a white 
dwarf plus low-mass companion, which results in a most 
likely spectral type of ^L6 for the companion. Based on the 
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Table 6. Summary of numbers at each stage of the process- 
ing. The columns are split for the spectroscopic and photometric 
methods. The spectroscopic sample is also further split by the 
classification of the optical spectra. 



Constraint 


Spectroscopic 


Photometric-only 




Objects 


Objects 


OetllolV lllti OUOO 






Colour Cuts 


7 AAA 




Spectroscopically Confirmed 






DA white dwarfs 


4636 




QSO 


1280 




1\TT UC 
IN 1 j I 1.3 


SAO 




Other white dwarfs 


661 




Other Objects 


27 




Objects cross matched with 






TTKTTm^ with detection in 

UlVli^OO W1L11 UCLCC L1U11 111 






any of V, J , H or K 


iyyu 


1 / I 1 


Spectroscopically Confirmed 






DA white dwarfs 


1275 




QSO 


328 




NLHS 


209 




Other white dwarfs 


172 




Other Objects 


6 




Objects with Detections in 






H &t K 


1075 


809 


K 


1108 


840 


Spectroscopically Confirmed 






DA white dwarfs 


571 




QSO 


316 




NLHS 


124 




Other white dwarfs 


94 




Other Objects 


3 




IR Excess Objects 


42 


105 



SDSS/UKIDSS data alone, it is impossible to distinguish 
between a low-mass companion and a dusty disc, but mid- 
infrared data data can break this degeneracy. We therefore 
classify infrared excess candidates that require a companion 
later than L7-8 as brown dwarf / dusty debris disc candi- 
dates 

SDSS J1228+1040 exhibits a very red H - K colour in 
the (z — H, H — K) colour-colour diagram shown in Fig. 1131 
where it is clearly separated from the white dwarf model 
sequence. This region of the colour-colour space is therefore 
likely to harbour white dwarfs with a strong 7f-band excess. 
A significant number of the IR excess candidates in Sect. 15.21 
also lie in this region, as well as somewhat below (corre- 
sponding toaJ? and Tf-band excess). SDSS J1228+1040 is 
also prominent in Fig.[]3]thanks to its relatively high proper 
motion and red H — K colour. In summary, it is encouraging 
that this benchmark system indeed stands out in the various 
diagnostics we have considered. 



6.1.2 SDSSJ1043+0855 

SDSS J1043+0855 is the second white dw arf in our sam- 
ple k nown to have a gaseous debris disc (Gansickc ct al. 
l2007f ). though the evidence for an infrared excess in the 



near- and mid-IR is marg inal (|Melis et alj|201fj Brinkworth 
et al. 2011, in prep.). Based on the published results on 
this object, we would not expect to detect any excess in the 
UKIDSS photometry. In our analysis of the SDSS spectrum, 
we find T cff = 17912±360 K and log(s) = 8.0 7±0.08, consis- 
tent w ithin the errors with the parameters in Gansickc ct al. 
(|2007h . The corresponding model is an excellent fit to both 
the SDSS and UKIDSS photometry, with no detection of 
an IR excess in any of the near-IR bands fFig. H6|l , Adopt- 
ing the K a flux from iMelis et al.l (|2010h rather than the 
UKIDSS measurement gives a 2a excess above our model, 
which would again not classify as an IR excess candidate 
within our classification scheme (Sect. [53]) . This conclusion 
is confirmed by the fact that SDSS J1043+0855 falls very 
close to the model white dwarf sequence in Fig. 1131 This 
system is an example where there either is no dusty disc, or 
where the dusty disc is too faint, e.g. edge-on, to be detected 
in the near-IR. 



6.1.3 SDSSJ1212+0136 

SDSS J1212+0136 is one of three magnetic (DAH) white 
dwarfs that were selected by our colour cut (Table[T} and 
that, when fitting their ugri photometry, show a substan- 
tial IR flux excess (Fig.fTTl Sect. 15.4.30 . This white dwarf 
w as first reported to have a magnetic field of ~ 13 MG 
bv lSchmidt et all (|2003T ) . ISchmidt et ail (|2005t ) subsequently 
detected a weak Ha emission line, from which they measured 
an orbital period of ~ 90 min. Based on t h e J b and magni- 
tude of SDSSJ1212+0136, ISchmidt et ail l|2005T ) concluded 
that the companion to the white dwarf is a brown dwarf. Ad- 
ditional studies in the near-IR confirmed the brown dwarf to 
have a spectral type in the range L5-L8, and detected vari- 
able cyclotron emission, indi cative of ongoing accretion onto 
the magnetic white dwa rf l|Debes et al.l 120061 ; iFarihi et al.l 
120081 ; iHowell et aLll2008T ). Observations at blue and ultra- 
violet wavelengths show a quasi-sinusoidal flux modulation 
interpreted as emission from an accretio n-heated polecap 
l|Burleigh et al.l l200fj| ; iLinnell et all l2010h which is typical 
of strongly magnetic cataclysmic variables ( Gansick e et al.l 
1 19951 ; lArauio-Betancor et al.l 120051 ; iGansicke et al.l l2006al l . 
Despite the fact that no state of high accretion activity 
has been observed in SDSS J1212+0136, all observational 
evidence suggests that is nearly a twin o f the prototypical 
magn etic cataclysmic variable EF Eri (e.g. iBeuermann et ah! 
I2000T ). 

Taken on its own, the very red H — K colour 
of SDSS J1212+0136 would suggests it to be a quasar 
(Sect. l5~4.2l Fig. Q3}. However, its stellar nature is confirmed 
by the detection of a significant and large proper motion 
(Fig. I14p . Recovering this DAH white dwarf with a brown 
dwarf companion shows that our methods to identify DA 
white dwarfs with IR-excess are sufficiently robust to also 
find non-DA white dwarfs with genuine IR excess. 



6.2 Example Spectroscopic IR Excess Candidates 

6.2.1 SDSS J0135+1U5 

SDSS J0135+1445 is a clear candidate for being a cool white 
dwarf with a late-type stellar companion or brown dwarf. 
An excess is seen to extend over all the UKIDSS bands in 
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Wavelength |>m] Wavelength [>m] 

Figu re 15. SDSS J1228+1 040 is one of two DA white dwarfs in our SDSS/UKIDSS sample that are known to hav e a gaseous debris 
disc llGansicke et al.ll2006bl) . It also exhibits a substantial infrared excess in VLT/ISAAC and Spitzer observations I IBrinkworth et al,l 
2009), and hence serves as a general benchmark for our method. Fitting the SDSS spectrum (left panel) results in T e ff = 22037 dz 199K, 
log(q) = 8.19±0.04, implying M w( j = 0.74±0.02Mq , and a distance of 134±3pc, consistent with the parameters derived by Gansickc ct al. 
(2006b). The SDSS spectrum is shown in gray, the best-fit model in black, observed ugrizY J H K fluxes in red, and model fluxes in these 
bands as black circles. The lower panel shows the residuals from the fit (in flux units on the left axis, and statistical significance on the 
right axis). The corresponding fit to the ugri photometry gives a somewhat lower temperature, T fj = 20000±|[qK, log(g) = 8.5±Q'JgK, 
which has no noticeable effect on the detection of the infrared excess (right panel). 




Wavelength |>m] Wavelength [/jm] 

Figure 16. SDSS J1043+0855 is the second of only two DA white dwa rfs in our SDSS/UK IDSS sample that is known to have a dusty 
debris disc. This ob ject was discovere d through optical Ca emission jGansicke et alj|2007t) and has at best a marginal flux excess in 
the Spitzer mid-IR jMelis et alj|20ld Brinkworth et al. 2011, in prep.). Fitting the SDSS spectrum results in T off = 17912 ± 360K, 
log(cj) = 8.07 ± 0.08, the corresponding fit to the ugriz photometry results in T eff = 16000±^gK, log(p) = 7.5±g;| 2 K (right panel). 
No IR excess is detected in the near-IR using either method, consistent with the inconspicuous location of SDSS J1043+0855 in the 
(z — H,H — K) and proper motion vs H — K diagrams (Figs . [T3l and Il4l respectively). 



Fig. [6] Modelling of the companion suggests that its spectral 
type is between L7 and L8. Spectral fitting implies a white 
dwarf of T cff = 7467 ± 18K at a distance of 69 ± 2pc. A 
T c ff = 8000±2(|K is calculated from fitting of the photome- 
try, broadly similar to the recently d iscovered (resolved ) DA 
plus brown dwarf binary PHL5038 jSteele et al.ll2009r i. 



13432 ±710K, log(g) = 7.81 ±0.15, with an implied distance 
d = 349±32pc and white dwarf mass M wd = O.5O±O.O8M . 
Fitting of the companion type was inconclusive, but the 
shape of the SED is similar to that of the benchmark object, 
SDSS J1228+1040. Therefore SDSS J0753+2447 is classified 
as a brown dwarf or disc candidate. 



6.2.2 SDSSJ0753+244 7 

SDSS J0753+2447 is a very strong candidate for being a DA 
white dwarf with a late-type brown dwarf or debris disc 
(Fig.[[8]). The fit to the SDSS spectrum implies a T cfi = 



6.2.3 SDSS J mi +1035 

SDSS J1247+1035 is a candidate for having a brown dwarf 
companion or dusty debris disk, however the UKIDSS K- 
band is only in excess by ~ 3cr over the white dwarf model. 
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Figure 17. SDSS J1212+0136, a short-period binary contain- 
ing a magnetic DA(H) w hite dwarf plus a brown dwarf compan- 
ion ([Schmidt et alj|2005h . The system is undergoing weak mass 
transfer, producing cycl otron emission that contributes to the 
observed near-IR excess llDebes et alj|2006t iBurleigh et alj|2006t 
iFarihi et afll2008h . SDSS J1212+0136 was picked up by our fit to 
all the photometric objects satisfying the colour cuts designed to 
find DA white dwarfs (Tablc[T). 
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Figure 18. SDSS J0753+2447, a DA white dwarf plus dusty 
disc or low-mass companion candidate. The best-fit to the SDSS 
spectrum (black line) gives T c g = 13432 ± 710K and log(p) = 
7.81 ± 0.15 at a distance of 349 ± 32pc. The implied mass of the 
white dwarf is 0.50 ± 0.08Mq. Fitting of the 4cr if-band excess 
with main sequence star models proved inconclusive, however, 
considering the similarity to SDSS 1228+1040 (Fig. \TE^, this is an 
excellent candidate for having a brown dwarf companion or debris 
disc. Fitting the ugri photometry leads to T off = 12000±^°K 
and log(g) = 7.75±q |j at 327}|pc, with no change to the con- 
clusions as to the nature of the excess. 

Far-IR photometry of the object is required to confirm the 
IR excess. 

6.2.4 SDSSJ1557+0916 

UKIDSS photometry of SDSS J1557+0916 shows a 4a K 
band excess for both the spectroscopic and photometric fit- 
ting methods. The spectroscopy and photometric T e fi differ 
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Figure 19. SDSS J0959-0200, a photometric-only DA white 
dwarf candidate (T cff = 120Q0±jjJ{j°K, log(ff) = 8.00±J;§S?) that 
exhibits a strong i<"-band excess, making it a strong candidate 
for having either a late-type brown dwarf companion or a dusty 
debris disc. 

by 3800K, however, this does not significantly affect the re- 
sult. This is a good example of where reddening is probably 
reducing the blue flux. Reddening the white dwarf model 
spectrum by E(B — V) ~ 0.05 brings the overall SED in 
line with the SDSS optical spectrum. Modelling of the com- 
panion object proved inconclusive as to its spectral type 
(Table©. SDSS J1557+0916 is a good candidate for having 
a dusty disc or low-mass companion based on its IR spectral 
shape. 

6.2.5 SDSS J 2220- 004 1 

PHL5038 ( = SDSS J2220-0041) is a wide (0.94") binary 
containing a cool (~ 8000 K) white dwarf with an ~L8 com- 
panion, only the fourth wh ite dwarf plus brown dwarf binary 
known (jSteele et alj|2009h . 

6.3 Example Photometric IR Excess Candidates 

As described in Sect. l5.ll2l we fitted all photometric objects 
satisfying our DA white dwarf constraint set (Table[l} with 
DA model spectra, independent on whether they also have 
an SDSS spectrum. 

6.3.1 DA white dwarf candidates 

6.3.1.1 SDSS J 0959 0200 SDSS J0959-0200 is a 
photometric-only DA candidate with T e g = 12000±mo° K 
and log(g) = 8.00±J:i. The UKIDSS Tv-band magnitude 
shows a large (5<r) excess over the white dwarf model, no 
excess is seen at shorter wavelengths. This object is a strong 
candidate for having a very late type brown dwarf compan- 
ion, or a dusty debris disc. 

6.3.1.2 SDSS J '1221+1245 A second interest- 
ing photometric-only DA white dwarf candidate is 
SDSS J1221+1245. The white dwarf is best fitted by a 
model with system parameters T e ff = 12000±27o° K and 
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Figure 20. SDSS J1221+1245, a photometric-only DA white 
dwarf candidate (T cff = 12000±^J°K, log(g) = 8.00+g;^) that 
exhibits a borderline i-f-band excess. It is good candidate for hav- 
ing either a late-type brown dwarf companion or a dusty debris 
disc. 



log( 5 ) = 8.00±o:|i The UKIDSS if -band magnitude shows 
a borderline excess over the white dwarf model, but no 
excess is seen at shorter wavelengths. This object is again 
a good candidate for having a very late type brown dwarf 
companion, or a dusty debris disc. 



6.3.2 Other Composite Objects 

Inspection of Simbad reveals that four of the photometric- 
only DA candidates with IR excess (DA:ire and DA:ire:, Ta- 
ble[8j) are previously known (pre-)white dwarf binaries and 
one pulsating subdwarf, which provides a preview on the 
mixture of objects that can be expected within this sam- 
ple. It also underlines that the method efficiently identifies 
genuine infrared-excess white dwarfs. 

6.3.2.1 SDSS J001 6+0704 =PG 0014+067 
iBrassard et~aT1 l|200ll) identified PG 0014+067 as a pulsat- 
ing sdB with Tea = 33550 ± 380K and log(g) = 5.77 ± 0.10. 
Fitting the ugri photometry with DA model atmospheres 
results in T cff = 24000±!oo°K with the surface gravity 
fixed to log(g) = 8, and reveals a clear IR excess. At such 
high temperatures, the slope of the optical and near-IR 
SED of this object is close to a Rayleigh- Jeans distribution, 
and while modelling the photometric data with DA models 
may not be perfect, we believe that PG 0014+067 does 
exhibit a genuine IR excess. I n their asteroseismological 
analysis, IBrassard et al.l |200lJ ) found that the pulsation 
frequency spectrum of PG 0014+067 exhibits fine structure 
that they tentatively interpre ted as a rotational period of 
29.2 ± 0.9 h, revised later by l|jefferv et al.ll2005l ) to ~ 4d. 
One possibility is that PG 0014+067 has a close low-mass 
binary companion with an orbital period of a few days and 
that the sdB is tidally locked, rotating at the same period. 
This hypothesis can be tested by a radial velocity study of 
this subdwarf. 



6.3.2.2 Cataclysmic Variables BKLyn 

( = SDSS J0920+3356) and HS 0139+0559 



f = SDSS0141+0614) are novalike v a riables 
jPobrzvcka fc Howelll Il992l ; lAungweroiwit et all I2005T I 
with optically thick accretion discs, and their optical 
colours are similar to that of hot white dwarfs or subdwarfs. 
However, their companion stars and cooler outer regions of 
the accretion discs start to contributed noticeably in the 
near-IR. 

6.3.2.3 Detached Binaries Abell31 

( = SDSSJ0854+0853, PN G219.1+31.2) is a planetary 
nebula with a nearby (0.26") M- dwarf, both stars are most 
likely an associated wide binary (|Ciardullo et al.lll999f ). 

GKVir ( = SDSS J1415+0117) is an eclipsing binary 
containing a hot (~ 48800 K) white dwarf pl us an ~M4V 
companion with an orbita l period of 0-344 d dGreen et al.l 
ll978i ; lFulbright et al.lll993l ; iParsons et alj|20ich . 



7 COMPARISON WITH SDSS DR6 WHITE 
DWARF-MAIN SEQUENCE BINARIES 

iRebassa-Mansergas et all |2010l . hereafter RM10) compiled 
a catalogue of white dwarf-main sequence (WDMS) bi- 
naries from all spectroscopic objects within SDSS DR6. 
Given that their detection method was based on optical 
data alone, RM10 were primarily sensitive to white dwarfs 
with M-type companions. The distribution of their WDMS 
binaries as a function of effective temperature and com- 
panion star spectral type shows a clear concentration of 
T c ff = 10 — 20kK white dwarfs with ~M4-type compan- 
ions (Fig. 1211 left panel). The large luminosity of hot white 
dwarfs prevents the identification of low-mass companions 
around them, explaining the relative dearth of late spec- 
tral types at higher temperatures. The decreasing number 
of very late M-dwarfs (>M6) could also be affected to some 
degree by the same contrast problem, however, it is known 
that the companion mass distribution of WDMS binaries is 
dropping towards the low-mass end of the main sequence 
(|Farihi et al.ll2005T l. 

The sample of WDMS binaries of RM10 provides a nat- 
ural comparison for the work done here. We have subjected 
their entire sample of WDMS binaries to our DA colour 
cuts (Table[T]), finding that only 93 of the 1602 systems fall 
within the colour cuts. This small number is not too sur- 
prising, as the red flux from the M-dwarf companions moves 
the majority of RMlO's WDMS binaries out of our colour 
selection. RM10 list WD temperature and companion spec- 
tral type for 53 of these 93 systems. In contrast to our work 
here (Sect. [5T2|) . the study of RM10 included (partially) re- 
solved systems. Consequently, we removed 21 that appeared 
resolved in the SDSS (or UKIDSS images where available), 
which leaves us with 32 objects in common. Figure [2T1 shows 
that the two samples only overlap for systems where the 
companion is relatively faint in the optical compared to the 
white dwarf, which is expected as our DA selection needs 
the white dwarf to dominate. Finally, 10 of the 32 objects 
are in the UKIDSS footprint. 

We cross-correlated the white dwarfs with an IR ex- 
cess from our spectroscopic and photometric samples (Ta- 
bleland©, and the WDMS catalogue of RM10 and found 
19 objects in common (Table [T0)l . This is comprised of the 
10 objects above that we expect to be in the sample, along 
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with 9 others where no spectral type is listed in RMlO's cat- 
alogue. There is in general a good agreement between the 
white dwarf system parameters, however the underestima- 
tion of photometric temperatures is highlighted again. The 
spectroscopic sample and the WDMS catalogue largely find 
companion types within two spectral types of each other. 
All of these are M-type companions as expected from the 
sensitivity of the WDMS catalogue. Five objects are found 
to have an excess in the photometric method, but were re- 
jected from the spectroscopic sample because signatures of 
a main sequence star companion can be seen in the opti- 
cal spectra, whereas the spectroscopic sample only contains 
pure DA white dwarfs. Those are marked as such in Table flOl 
and suffer from the same companion type biases discussed 
in Section l5"XTl 



8 CONFIRMATION OF IR EXCESS 
CANDIDATES IN WISE 

We cross correlated all spectroscopic and photometric-only 
IR excess candidates (from Table[7] and [8] respectively) with 
WISE PDR within 2.5". Of the 42 IR excess objects from 
the spectroscopic method, 3 had a detection in at least one 
WISE band. Similarly 14 of the 67 photometric-only objects 
were detected and are listed in the online master table. 

The three spectroscopically confirmed DA white dwarfs 
with IR excesses and WISE data: SDSS J0236-0103, 
SDSSJ0847+2831 and SDSS J1448+0713, are all confirmed 
to have an excess in the WISE 3.4 and 4.6/im bands. They 
are however all predicted to have late M-type companions 
and are not brown dwarf or debris disk candidates. 

Of the 14 photometric-only IR excess candidates with 
WISE data, SDSS J1524-0128 and SDSS J1549+0325 are 
not found to have an excess in the 3.4 and 4.6pm bands. 
Similarly, because we do not trust the effective tempera- 
ture of the white dwarf fit, we also do not believe the far- 
IR excess found for SDSS J0841+0501, SDSS J1441+0137, 
SDSSJ1538+0644 and SDSS J1551-0118. This is indicated 
by a flat, constant excess over the white dwarf model and 
is continued into the far-IR WISE data in these cases (see 
Fig.|22]|. 

The remaining 8 objects (SDSS J0207+0715, 
SDSSJ0742+2857, SDSS J0751+2002, SDSS J0920+3356, 
SDSSJ1448+0812, SDSS J1456+1040, SDSS J1538+2957, 
and SDSS J1635+2912) are confirmed to have an IR excess 
consistent with a late type companion in the WISE far-IR 
data. Some interesting examples of these are shown in 
Fig.rjS] SDSSJ1538+2957 (Fig.[23J is predicted to have 
an M8-type companion from the photometric method. 
However, the spectral shape of the excess is found to be 
inconsistent with such an early type companion. The excess 
in UKIDSS, and now WISE, is more indicative of a later 
type brown dwarf companion or dusty debris disk. This 
mismatch is most likely caused by over estimating the 
distance to the white dwarf in the photometric method. 
SDSSJ1635+2912 (Fig.123)) is one of the photometric-only 
debris disk candidate systems. The WISE 3.5pm flux 
confirms the infrared excess, but still leaves the origin of 
the excess, brown dwarf or disk, open. 




Wavelength (um) 



Figure 22. SED of SDSS J1538+0644. The SDSS, UKIDSS and 
WISE 3.4 and 4.6/nn fluxes are shown as red circles. The best 
fit white dwarf model to the SDSS photometry is shown as a 
black line. The shape of the IR excess is not consistent with any 
companion or disk. It is more likely that the excess is an artifact 
of overestimating the white dwarf effective temperature. 

9 DISCUSSION 

9.1 Selection of DA white dwarfs from optical 
photometry 

We have developed dedicated colour-colour cuts to select 
DA white dwarfs only from their optical photometry (Ta- 
bled Fig. [2]). This method can easily be optimised to priori- 
tise for completeness, efficiency or a compromise of both. We 
have demonstrated that a high completeness (95.4%) can be 
obtained with a reasonable efficiency (62.3%) based on the 
ugriz data from SDSSDR7. The strengths of this approach, 
is that it provides substantially larger and statistically bet- 
ter characterised white dwarf samples (Scct. [3~3l when com- 
pared to sp ectroscopic catal o gues such as iMcCook fc Sionl 
(1999) and lEisenstein et al.l (|2006t ). We have also investi- 
gated methods to account for the contamination of the DA 
white dwarf candidate sample by quasars and NLHS. This 
method can be adapted to obtain large uniform samples 
of white dwarfs from other multicol our optical photo met- 
ric surveys, such as e.g. SkyMapper <|Keller et alj|2007ft . 

9.2 Implications for the brown dwarf desert and 
the numbers of white dwarfs with dusty 
debris discs 

Fitting the optical spectroscopy or photometry and probing 
for IR flux ex cess above the best-fit white dwarf, as previ- 
ously done bv lTremblav fc Bergeron! (120071 ). has proven to 
be an efficient and robust approach. It also allows one to 
identify hot white dwarfs with low-level excesses which have 
IR colours that are very similar to the bulk population of 
white dwarfs (Fig. ll3[l . Our search is sensitive to companions 
as late as ~L8, and to warm dusty debris discs. 

Of the 1275 spectroscopically confirmed DA white 
dwarfs with at least one UKIDSS magnitude, 26 (2.0%) are 
found (or are candidates) to have an IR excess consistent 
with a LO-type companion or later, i.e. a brown dwarf. The 
exact cut off of where the brown dwarf sequence starts is age 
dependent and can vary from mid-M to mid-L-type. Taking 
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Figure 21. The distribution of the WDMS binaries from the catalogue of Rcbassa-Manscrgas ct al. (2010) as a function of effective 
temperature of the white dwarf and spectral type of the companion star. The left hand panel shows all the 1173 objects with values for 
both T e ff and log(g). 53 of these are contained within the colour-colour region defined in Table[l] Excluding spatially resolved binaries 
leaves 32 objects satisfying all the criteria of our DA white dwarf selection, which are shown in the right hand panel. 
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Figure 23. SEDs of SDSS J1538+2957 (left) and SDSS J1635+2912 (right). The plot follows the same format as Fig.[22] SDSS J1538+2957 
and SDSS J1635+2912 are predicted herein to have an M7: and an L6:-type companion or potential debris disk respectively. 



the white dwarf mass distribution and the initial-final mass 
relation we can calculate the white dwarf progenitor masses. 
This with the white dwarf cooling age allows us to estimate 
the total age of the system. The average total age of the com- 
panion is a few Gyr, and in this regime LO-type is a suitable 
cut-off for being a brown dwarf. Taking only the systems 
where we are confident of the excess (not "DAire:" in Ta- 
ble[7|) gives a lower limit for the number of white dwarfs with 
brown dwarf type companions of 0.8%. This is compatible 
with previous estimates. A n adaptive optics imaging sur - 
vey of 266 solar- like stars bv lMetchev fe Hillenbrand! j2009l ) 
found a w ide (28-1590 AU) su b-stellar companion fraction of 
3.2 + ^%. lFarihi et al l (2005) carried out an extensive near- 
IR imaging survey for both wide and unresolved low-mass 
companions to 394 known white dwarfs, finding an overall 
stellar companion fraction of 22%, and a brown dwarf com- 
panion fraction of < 0.5%. 

When using the photometric fitting method, 105 of 
the 1771 photometric-only DA white dwarf candidates with 
UKIDSS data exhibit an IR excess, of which we eliminate 38 



likely quasars (Sect. [57131 and see online QSO table). Taking 
into account our estimate for the contamination by NLHS 
(Sect. l5"T4"73")l . we find that ~ 2.7% of white dwarfs in the 
photometric-only sample have an IR excess, of which 1.8% 
are candidates for having a brown dwarf companion. This 
is consistent with the frequencies found in the spectroscopic 
sample, though not as secure. 

Taking only the objects with a detection in the K-b&nd 
(required for detecting a disk), of the 571 spectroscopically 
confirmed DA white dwarfs, 7 (1.2%) are found to have an 
IR excess compatible with a companion spectral type later 
or equal to L8, and are hence viable disc candidates. Simi- 
larly, 0.5% (4/840) of the photometric-only sample are disc 
candidates, or, scaling for the 62.3% efficiency of our DA 
white dwarf s election, 0.8% (4/52 3). This is consistent with 
the results of iFarihi et alj (120091 ), who estimated that the 
frequency of white dwarfs with dusty debris discs is at least 
1%. Any search for debris discs based on JC-band data is 
biased towards warm, bright circumstellar dust rather than 
faint discs or narrow rings. Of the known white dwarfs with 
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a dusty debris disc, only about half reveal themsel ves short- 
ward of 3 urn (|Kilic et al.ll2006l ; iFarihi et al1l2009l ). 

To confirm the IR excess candidates we have cross 
matched all the spectroscopic (Table[7)l and photometric- 
only (Table© candidates with the far-IR WISE PDR. This 
has provided 3.4 and 4.6/zm fluxes for 7% (3/42) and 21% 
(14/64) of the spectroscopic and photometric-only IR excess 
candidates respectively. We find that all three of the spec- 
troscopic IR excess candidates also have an excess in the 
far-IR. 

However, in the photometric-only sample, a total of six 
objects are found not to have a definite IR excess at > 3^tm. 
The remaining 8 of 14 white dwarf candidates have a real 
IR excess in the WISE fluxes. We therefore find that ~ 60% 
of the photometric-only IR excess candidates have real IR 
excesses consistent with a MS star companion, brown dwarf 
companion or dusty debris disk. 

Determining the true nature of the white dwarfs with an 
excess in K only will require far-IR data to distinguish be- 
tween a late-type brown dwarf companion or a dusty debris 
disc. More generally, deeper near-IR and/or far-IR follow 
up observations will be required to verify, or refute, the IR 
excess candidates (DAire: and DA:ire:) , as some of these can- 
didates are likely to be spurious. Similarly, while our colour 
selection (TableQ]) has a high efficiency, optical spectroscopy 
of the white dwarf candidates (DA:ire and DA:ire:) will be 
necessary to confirm their (DA) white dwarf nature. 



10 CONCLUSIONS 

We have developed a detailed method to select DA white 
dwarfs in ugriz colour space, and applying this method to 
SDSSDR7, we have identified 7444 DA white dwarfs with 
g ^ 19 and SDSS spectroscopy, approxi mately 70% more 
than the corresponding number from DR4 (|Eisenstein et all 
2006), and 9341 photometric-only DA candidates. Using this 
sample, we estimate the spectroscopic completeness of DA 
white dwarfs with T > 8000 K in SDSSDR7 to be ~ 44%. 

We cross-correlated these spectroscopic and photomet- 
ric DA samples with UKIDSS DR8 to carry out the currently 
largest and deepest untargeted search for low-mass compan- 
ions to and debris discs around DA white dwarfs. This search 
led to the identification of a significant number of DA white 
dwarfs with low-mass companions, including several brown 
dwarf and dusty debris disc candidates. Similar studies to 
the one presented here making use of the full UKIDSS area, 
as well as the corresponding surveys in the southern hemi- 
sphere, e.g. Skymapper, VISTA and VHS, promise to further 
increase the sample sizes of these types of systems. 



NOTE 

Before submission the authors became aware of the study 
of Steele et al. (2011, in prep.) who have carried out 
an independent study, cross-m atching the spect r oscop ically 
confirmed white dwar fs from iMcCook &: Sionl (1 19991 ) and 
lEisenstein et al.1 (|2006T ) with UKIDSS DR8. Comparison of 
their results with our spectroscopic sample shows that both 
studies lead to broadly consistent results. 
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Table 7. 42 DA white dwarf IR excess candidates from the spectroscopic method. In all columns, ':' indicates an uncertain classification. 
The notes column indicates the level of certainty of the IR excess. The key to the notes2 column is: bg (background contamination), 
disc (disc candidate), moreIR (more IR data required ), resolved (resolved or partially resolved binary). The "Eisenstein" column is the 
classifications according to the Eisenstein ct al. (2006) catalogue. The column "JHK" firstly indicates whether the object has a UKIDSS 
J, H and/or K magnitude. Secondly, if the corresponding letter is bold, it means an excess was found in that band. 



Coord T cff (X) log(ff) M wd (M ) R wd (10 8 cm) d(pc) Companion Notes Notes2 Eisenstein JHK 



00324-0739 


21 045 + 240 


7 42 4- o 05 

I .rlL _l_ U.UU 


37 4- o 02 

VJ.iJi _1_ U.Uii 


1 37 4- 05 

1 . U I _l_ u . uu 


308 + 1 5 

JC/U _l_ 1U 


M6 


DAire: 


1 LoUlVLU, '-^o' 


DA 


JHK 


0039 — 0030 

V ' V ' U >J V J U U U 


12392 ± 544 


7 35 + 1 6 

1 .UU _1_ 


32 + 05 


1 37 + 1 5 

1 . U 1 _1_ U . 1 u 


417 ± 44 


M9 


DAire: 




DA antn 


JHK 


01 354-1 445 

U 1 JO | IrtrlU 


7467 ± 18 


7 34 4. o 04 

I . Urt _|_ u.Urt 


29 + 01 

U.iC/ _1_ U.U1 


1 33 4- 03 

l.uu _l_ U.UU 


69 + 2 


L6 


DAire 




DA 


JHK 


0207+0702 

Uiiu i n^u i ul 


1 0073 + 77 

1UVJ 1 U _l_ i ( 


a 32 + 00 

O.UL _l_ U.U3 


80 + 06 

U.UU _1_ U.UU 


72 + 05 

U. 1 L _l_ U.UU 


162 ± 10 


L7 


DAire: 






JHK 


n2Sfi oi 03 

UiJU U 1 U U 


20566 + 4Q8 


7 fit; 4. o 00 


4fi 4- o 04 

U.rtU _1_ U.Url 


1 1 5 4. o 07 

1 . 1U _l_ U.Ul 


411 ± 25 


M6 


DAire 




DA+M: 


JHK 


0253—0027 


1 8Q74 + 286 

lOi? 1 t: _l_ LOU 


7 71 4- o 06 

i . i i _i_ u. uu 


47 4- o 03 

U.rt I _1_ U.UtJ 


1 114.005 

1.11 _l_ U.UU 


300 + 1 5 

UOU _l_ 1U 


L i 


DAire 


Tnr^fi^TR 

111U1 L11V 


DA antn 

1 J i \ < 1 ULU 


JH 


07484-2058 

U 1 rtOn^LUUO 


8fi726 4- 7732 

OU I LU _l_ 1 I U L 


7 1 5 4. o 27 

1 . 1U _l_ U. L J 


53 4- o 06 

U.JJ _1_ U.UU 


2 23 + 1 37 

. LJ _l_ l.U 1 


2021 + 469 

LUL1 _1_ lUJ 


M3 


DAire 


1 CdWI v cu • 


DA+M: 


JHK 


07534-2447 

U I JOT"" I 


1 3432 + 71 

1 JUii _1_ 1 1U 


7 »1 4-015 

I . O 1 _l_ U . 1U 


50 + 08 

U.UU _1_ u.uo 


1 02 + 1 

l.vij _l_ U. 1U 


349 ± 32 


L5: 


DAire 




DA oiifn 

TJ T± — Ch Ll LU 


JHK 


08474-2831 

UOrt I |^ iOUl 


1 2828 + 030 


7 7« 4. o 20 

I . i u _i_ u. lu 


48 + 1 1 

U.rtO _1_ U.ll 


1 05 + 014 

l.UU _l_ U.lrfc 


357 ± 42 


M6 


DAire 




DA+M: 


JHK 


09074-0536 

KJUKJ 1 1 UUUU 


1 041 fi + 458 

li/'rllU _l_ tuU 


7 7a 4. o 00 


n 51 4- o 05 

U.U1 _1_ U.UU 


1 06 + 06 

l.UU _l_ U.UU 


367 ± 21 


L6 


DAire: 


Tin /"irpTR 
111U1 L^ll v. 


DA aiitn 

lvA_aU LU 


JHK 


09334-3200 


1 1 5K5 4- 381 


a 3« 4. o 20 

O.UU _l_ U.LU 


83 + 1 3 

U.OtJ _1_ U. ItJ 


70 + 1 1 

U. 1 U _l_ U.ll 


229 ± 34 


Li 


DAire 


1 toWl v • 




JHK 


09504-01 1 5 

U U U U I U 1 1U 


21 785 + 3fi5 

L 1 1 OU _l_ UUU 


7 on 4- 06 

1 . O tJ _|_ U . 


57 4- o 03 

U.Ul _1_ U.UU 


qo 4- 04 


329 ± 14 


L8 


DAire: 


111U1 U>11\, 


DA 


JH 


1 0024-0939 


21 785 + 808 

L 1 i OU _l_ OUO 


7 Q2 + 1 4 

1 . C/L _l_ (J. Irt 


58 + 08 

U.UO _1_ u.uo 


07 4- o 00 


515 ± 45 


L0 


DAire 




DA ?mtn 

J ' - \ < 1 ULU 


JHK 


1 01 04-0407 


1 3588 + 668 

1UUOO _l_ uuo 


7.76 ± 0.11 


48 + 05 

U.rtO _1_ U.UU 


1 05 + 07 

1 . UU _l_ U.Ul 


280 + 1 8 

LOU _l_ lO 


L8: 


DAire: 


111U1 U>11V 


DA aiitn 

J ' - \ < 1 ULU 


JH 


1 01 54-0425 


34526 + 86 


7 38 + 07 


n 41 4- 02 

V 7 . i 1 1 U.UL 


1 51 + 08 

1.U1 _1_ U.UU 


691 + 39 


L4 


DAire: 


lllOrellr* 


DA_auto 


JHK 


1 01 64-0020 

1U 1 U 1 UULU 


21 045 + 703 

L 1 ut:J _l_ 1 UU 


8 48 + fl 12 
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U.fiJ _1_ U.Ul 


fi4 4- o 06 

U.Urt _l_ U.UU 


325 4- 32 

ULU _l_ UL 
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1UU I 1 Ulu J 
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fio + 07 

U.UC7 _l_ U.Ul 


151 ± 14 
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1 1 204-0630 

J- 1 L U |^ UUtJCI 


20331 + fiOO 

iUOJl _l_ UU -J 


7 «0 + 1 1 

I .U3 _l_ U.ll 


47 4- o 05 

U.rt I _1_ U.UU 


1 1 3 4. o 00 

1.1U _l_ u.uc? 


517 ± 38 
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DAire: 


rPGi^l"vrf>H no" 

1 LoUlVLU, 


DA aiitn 

J ' - \ < X ULU 


j 


1 1 41 4-0420 


1 1 235 J- ai a 

llOUU _l_ OlO 


7 37 4. o 22 


32 + 08 

U. Ji _1_ U.UO 


1 35 4. o 21 

l.UU _l_ U . 1 


317 ± 46 


M7 


DAire 




DA+M: 


JHK 


1 2084-061 


2^076 + 1125 

ii JU 1 U _l_ lli J 


7.64 ± 0.17 


4fi 4- o 08 

U.rtU _1_ U.UO 


1.18 ± 0.14 


725 ± 81 


M7 


DAire 




DA aiitn 
TJJrlL_cb U LU 


JH 


1 21 84-0042 


1 1 1 73 + 1 SO 

111 1 O _l_ lOU 


c 30 4- o 1 1 

O . <J \J _l_ U.ll 


70 4- o 07 

U. 1 C7 _1_ U.Ul 


73 4- o 06 

U. 1 U _l_ U.UU 


171 ± 13 


L8 


DAire: 


i i inrr 1 R 

111U1 Lll v 


DA 


JH 


1228+1040 


22037 + 199 

— lu* y j _i_ t *j <J 


8 1 9 + 04 


74 + 02 

V 7 . 1 i 1 U.UL 


80 + 02 

U.OU _1_ U.UL 


134 ± 3 


L6: 


DAire 




DA antn 

\ y - \ II ULU 


JHK 


1 2464-0707 


1 0793 + 1 89 

1U 1 <?U _l_ 103 


c 03 + 1 6 


fi2 + 1 

U.UiJ _1_ U.1U 


88 + 1 

U.OO _l_ U. 1U 


200 + 21 

LUU 1 LI 


L2 


DAire 






HK 


1 2474-1 035 

J. L*± I 1UUU 


1 vol 2 + 1 50 


7 «2 4- o 04 

i .Oil _l_ vj-yj'-t 


52 + 02 

U.Uij _1_ U.UiJ 


1 03 + 03 

l.UU _l_ U.UU 


394 ± 9 


L6: 


DAire: 


rio" - roQ/^ilArPri ■ 

Uc, ,1 LoUl VLU. 




JHK 


13144-0057 


1 7707 + 1 6,3 

ii i u i _1_ iv J 


7 80 + 04 


,51 + 02 

U.Ul 1 U.UiJ 


1 04 + 03 

1 . U^r _1_ U.UU 


312 ± 8 


L6 


DAire: 






HK 


1 3204-001 8 

1 U L u \^ UUlu 


1 01 03 + 213 

ISlJJ _l_ ilj 


a 40 4- 04 


87 + 03 

U.Ol _1_ U.UU 


R8 + 02 

U.UO _l_ U.UL 


142 ± 5 


L8 


DAire: 


disk- 


DA 


JHK 




i qaqo J- 979 

lurlUL 7JZ L ( L 




n 4. n 09 


± . ul u.uu 


209 + 8 


M7 


DAire 




DA mitn 
JJxi_ctU LU 


JHK 


i QQi 4-0040 

1UU 1 1 UUt:U 


1 5Qf;4 4. 71 4 

IUC/Ut: _1_ 1 Irfc 


7.72 ± 0.17 


47 4- o OS 

U.rt I _1_ U.UO 


i 0Q + 1 2 

1.UC7 _l_ U.1L 


4«5 4. aq 

rfcUU _l_ rtO 


? 


DAire: 


Tin r\vc\Vi 

111U1 L^ll V 


DA 




1 341 4-0056 

IUtI |^ KJ U U U 


1 8330 + 436 

lOtJOU _l_ rtUU 


7 QQ _i_ o 1 o 

1 • tJa _l_ U . 1U 


fii 4- o 06 

U.Ul _1_ U.UU 


no + 06 

U . C/L _l_ U.UU 


3«0 + 22 

UUU _l_ LL 


M7 


DAire 






JHK 


1352+0910 


36154 ± 722 


7.49 ± 0.11 


0.45 ± 0.04 


1.39 ± 0.13 


849 ± 76 


M5 


DAire 






JHK 


1410+0225 


11565 ± 246 


8.63 ± 0.13 


1.00 + 0.07 


0.56 ± 0.06 


152 ± 16 


3= L8 


DAire: 


moreIR 


DA 


JH 


1425-0013 


10670 ± 169 


7.97 + 0.17 


0.58 + 0.10 


0.91 ± 0.10 


242 ± 26 


L8: 


DAire: 


moreIR 


DA 


JHK 


1448+0240 


14728 ± 247 


7.41 ± 0.06 


0.34 + 0.02 


1.34 + 0.06 


292 ± 13 


Li 


DAire 


moreIR 


DA 


JH 


1448+0713 


12250 ± 256 


7.87 + 0.09 


0.53 + 0.05 


0.98 ± 0.06 


143 + 8 


M8 


DAire 






JHK 


1450+0000 


17106 ± 175 


7.93 ± 0.05 


0.58 + 0.03 


0.95 ± 0.03 


311 + 9 


L8: 


DAire: 


moreIR 




JHK 


1557+0916 


22811 ± 734 


7.75 ± 0.11 


0.50 + 0.05 


1.09 ± 0.08 


488 ± 35 


LI: 


DAire 




DA_auto 


JHK 


1609+2905 


8773 ± 91 


8.26 ± 0.14 


0.76 ± 0.09 


0.75 ± 0.08 


122 ± 12 


^ L8 


DAire 






JHK 


1619+2533 


25595 ± 1139 


7.21 ± 0.15 


0.33 + 0.04 


1.65 ± 0.20 


879 ± 101 


M5 


DAire 






JHK 


1625+3026 


72136 ± 4550 


7.77 + 0.19 


0.64 + 0.08 


1.21 ± 0.21 


739 ± 112 


M5 


DAire 




DA_auto 


JHK 


2220-0041 


7467 ±21 


7.66 ± 0.12 


0.41 + 0.06 


1.10 + 0.08 


69 + 5 


^ L8 


DAire: 


resolved 


DA_auto 


JHK 


2225+0016 


11045 ± 235 


8.44 + 0.16 


0.88 + 0.10 


0.65 ± 0.08 


167 ± 20 


L6 


DAire: 


resolved: 


DA_auto 


JHK 


2331+1342 


11173 ± 284 


8.79 ± 0.15 


1.09 + 0.08 


0.49 ± 0.07 


141 ± 19 


L6 


DAire 




DA+M: 


JHK 



Table 8. 67 photometric-only IR excess candidates which are found to have an IR excess in the photometric method. The columns follow 
the same format as that in Table[7] The key to the "Notes2" column is: bf (visually bad fit), bg (background contamination), disc (disc 
candidate), hot WD (high T c ff white dwarf and so small changes in, e.g. flux calibration, remove the excess), moreIR (more IR data required), 
resolved (resolved or partially resolved binary), qso (quasar), ? (unknown). Also, the "Simbad Name" and "Simbad Class" columns give an 
indication of what is already known about the objects. If the surface gravity could not be constrained, it was assigned the canonical value 
of 8.0 and therefore no error could be defined. We do not quote the distances to the systems because the uncertainty is so large. The two 
objects marked as "QSO:" are borderline cases of the objects found in the online QSO table and are therefore left as candidates for being 
DA white dwarfs with an IR excess. 



Coord 






Companion 


Notes 


Notes2 


JHK 


Simbad Name 


Simbad Class 


0016+0704 


24000± 2 i!i° 


8.00 


L0 


DA: ire 




JHK 


V* EK Psc 


Variable Star 


01044-1459 


12000+i 1 , 


7 75+S JI 


L5 


DA: ire 


moreIR 


JH 






0141+0614 


nooo±TS,„ 


8.50+2/ PJ 


L8: 


DA: ire: 


bf 


JHK 


2MASS 


Nova 


0207+0715 


14000±?i 4 P, 

1940 


9 25+HZ 


L6: 


DA: ire 


bf 


JHK 


PB6668 


Star 


0742+2857 


9000+™ 


7 25+2'i 2 . 


L6 


DA: ire 




JHK 






07514-2002 


20000+ 24e jS 


8 25+?-2§ 


L5 


DA: ire 




JHK 






0758+2406 


13000±|o° 


8 00±2 : ??, 


L2: 


DA: ire 




JHK 






0841+0501 


8000+JB2 


8.00 


^L8: 


DA: ire: 


bf 


JHK 






0842+0004 


22000+ i'-ii! 


8.00 


M9 


DA: ire 




JH 






08544-0853 


60000+S 890 


8.00 


M6 


DA: ire 




JHK 


PN_A66_31 


PN 


0906+0001 


7000±q° 


8.00 


^L8: 


DA: ire: 


bf 


JHK 






0920+3356 


12000±5*5o 


9 00±S ?l 


L8: 


DA:irc: 


bf 


JHK 


V* BK Lyn 


Nova 


0923+0652 


7000±o° 


8.00 


JsL8: 


DA: ire: 


bf be 


JHK 






0925—01 40 


1 4000+ 80 


7 50+°' 28 
'• OU=C 0.37 


M9 


DA:irc 




JHK 






09524-1205 


26000+i5 40 


8.00 


^L8: 


DA: ire: 


bf 


HK 






0959—0200 


12000+ 1160 


8 oo±H2 


L6 


DA: ire 




JHK 






0959+0330 


32000+22J!! 

o^ooozc 2 57o 


8.00 


L0 


DA:irc 




JHK 


PG0957+037 


UV 


1005—0225 


14000+ 5 2„ 


7 754-O.45 


L5: 


DA:ire: 


bf 


JHK 






10064-0032 


20000+ 4 '" 


9 ,50+0.91 


L8 


DA: ire 


bf: 


JHK 


PCI 004+008 


U V 


10384-1100 


10000+iZ° 


7 75±P iS 


L3 


DA: ire 




JHK 






1057+0628 


7000±f q 


8.00 


L8: 


DA: ire: 


QSO: 


JHK 






1 1034-0101 


7000+, 10 


8.00 


^L8: 


DA: ire: 


bf: 


JHK 


1VV20061 11 1 0336 34-01 01 41 




1103+1100 


34000± 2 .2?„ 


8.00 


L4 


DA: ire: 


hot WD, moreIR 


J 


PG1101 + 113 


SD 


1116+0755 


22000+f? 40 


8.00 


M9 


DA: ire 




JHK 






1153+0048 


15000+iiso 


9.50+0-36 


^L8: 


DA:ire: 


bf 


JHK 


LBQS1151+0104 


Star 


1221+1245 


12000±|^ 


8.00±0;5l 


L8: 


DA:ire: 


disk: 


JHK 






1224+0907 


7000±|3 


8.00 


L7 


DA:ire: 


bg:, resolved:, bf: 


JHK 






1237-0151 


20000± 4 2 1 


8.75±0;3| 


L4: 


DA:ire: 


resolved, bg: 


JHK 






1246+1600 


19000±|Jg 


9.50+ - 40 


^L8: 


DA:ire: 


bf: 


JHK 






1315+0245 


28000+jfgg 


8.00 


LI: 


DA:ire: 


resolved, bg: 


JHK 






1319+0152 


16000+^° 


8.25±g;|| 


>L8: 


DA:ire: 


bf 


HK 


WD1317+021 


DA 


1323+2615 


17000+1^° 


9.50+ - 89 


^L8: 


DA:ire: 


bf 


HK 






1334+0534 


30000±§|i° 


8.00 


^L8: 


DA:ire: 


bf: 


HK 






1348+1100 


22000±??° 


8.00 


LI 


DA: ire 


moreIR 


J 






1355+1454 


22000+22M 


8.00 


LI 


DA:ire 




JHK 


PB4150 


SD 


1416+1352 


38000±|°|°° 


8.00 


M6 


DA:ire 




JHK 


PK.003+66.1 


Galaxy 


1423-0138 


16000±|g° 


S-OO+oji 


^L8: 


DA:ire: 


moreIR 


JH 






1441+0137 


26000±|2, 


8.00 


L6: 


DA: ire: 


bf: 


JHK 


PG1438+018 


SD 



Table 8 - continued 



to 



Coord 


T eff (K) 


log(s) 


Companion 


Notes 


Notes2 


JHK 


Simbad Name 


Simbad Class 


1442+0910 


20000+JTQ 


7 75+1. J3 


M7 


DA: ire 




JHK 






1448+0812 


24000+ 1150 
^ uuu=c 2380 


8.00 


M7 


DA: ire 




JHK 


WD1 446+028 

VV 1 J 1 ^^y ) |W — (; 


DA 


1455+0458 


30000±||g° 


8.00 


M8 


DA: ire 




JHK 






1456+1040 


14000±} 4 §S 


8 75±S-S1 


LI: 


DA:ire 




JHK 






1 507+0724 


26000+2 4 98 


8.00 


M7 


DA: ire 




JHK 






1523+0023 


24000±1 3 ,§8 


8.00 


M8 


DA: ire 




JHK 






1524—0128 


18000±|S 


9 25±8?2 


^L8: 


DA: ire: 


bf 


JHK 


BPS_CS_22890— 0079 


SD 


1527+0802 


20000±i§ 3 ) ° 


9.50± - 63 


>L8: 


DA:ire: 




JHK 






1538+0644 


10000±|8 


7 25±2-i9 


L6 


DA:ire: 


bf: 


JHK 






1538+2957 


20000±f l 2 ' 


5 75±S'I5 

u. < .84 


M7: 


DA:ire 




JHK 






1540+2908 


34000± 4 2IS 


8.00 


M7 


DA: ire 
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8.00 


LI: 


DA:ire 


resolved: 


JHK 


2MASS 


Blue 


2308+0658 


28000±^6g 


8.00 


M6 


DA:irc 




JHK 






2326+1230 


18000±2°o 


Q 7C+0.58 


>L8: 


DA:ire: 


bf: 


JHK 


2MASS 


Star 


2330+1450 


13000±«g° 


9.50± - 64 


>L8: 


DA:ire: 


bf: 


JHK 






2335+1047 


11000±f§° 


7 75+O.55 
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Table 9. IR excess candidates which are found to have an excess in either the spectroscopic or photometric-only methods, but all have 
an SDSS spectrum for classification. The table follows the same format as Table[7] and [8] with addition that S and P in the titles refers 
to the spectroscopic and photometric methods respectively. Fiber and fit parameters are only shown if an IR excess was found in the 
corresponding method. 
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Table 9 - continued 
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1314+0057 


17707 ± 163 


7.80 ± 0.04 
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9 no+ 123 




>L8: 
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8.00+ 
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7467 ±21 
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8.00+ 
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DAire: 
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L6 


L5 
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DA 
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DA:ire: 
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15000±}|30 
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L8 
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Table 10. The 19 IR excess candidates from both the spectroscopic and photometric fitting methods that are also found to be white dwarf 
+ main sequence binaries in Rcb assa-Ma nscrgas ct al. (2010). Objects from the spectroscopic and photometric IR excess candidates are 
denoted as "S" and "P" in the "Method" column respectively. The first set of system parameters are from the methods described herein. 
Where both spe ctroscopic and photometri c IR excesses were found, the quoted parameters are those from the spectroscopic method, 
the second from Rebassa-Manscrgas et al.. In some cases, such as SDSS J0207+0702, the two spectroscopic fitting techniques disagree 
significant ly, this is where a differen t hot/cold fit solution was chosen. Photometric temperatures may also disagree significantly from 
thelRebassa-Manscrg as et al.l d201Cf) temperatures as dis cussed in Section[5TL2] The columns labelled with a "R-M" prefix are the best 
fit taken from Table 7 of iRebassa-Mansergas et all j201Ch . R-MSp=-l indicates that no spectral type could be associated. 



Name T cft log(ff) Method Companion Notes R-M T eff R-M log(g) R-M Sp 



0032+0739 


21045 ± 249 


7.42 ± 0.05 


S 


M6 


DAire: 


21045 ± 371 


7.43 ± 0.07 


M2 


0039-0030 


12392 ± 544 


7.35 ± 0.16 


SP 


M9 


DAire: 


12392 ± 737 


7.39 + 0.22 


-1 


0207+0702 


10073 ± 77 


8.32 ± 0.09 


SP 


L7 


DAire: 


18756 ± 502 


7.42 + 0.10 


-1 


0236-0103 


20566 ± 498 


7.66 ± 0.09 


SP 


M6 


DAire 


21289 ± 723 


7.83 + 0.13 


MS 


0748+2058 


86726 ± 7732 


7.15 ± 0.27 


SP 


M3 


DAire 


87731 ± 10728 


7.17 + 0.38 


M2 


0847+2831 


12828 ± 930 


7.76 ± 0.20 


SP 


M6 


DAire 


13588 ± 923 


7.76 + 0.18 


M6 


0851+0330 


20000+15* 


9.50+000 


P 


L6 


DA:ire 


22550 ± 520 


7.40 + 0.08 


M6 


1037+0139 


11433 ± 208 


8.38 + 0.13 


SP 


L5: 


DAire 


16717 ± 602 


7.78 + 0.14 


M5 


1106+0737 


28000±|060 


8.00 


p 


M5 


DA:ire 


36154 ± 540 


7.74 + 0.09 


M5 


1141+0420 


11835 + 818 


7.37 + 0.22 


SP 


M7 


DAire 


12110 + 780 


7.39 + 0.26 


-1 


1329+1230 


13432 ± 272 


7.42 ± 0.05 


SP 


M7 


DAire 


12250 ± 1036 


7.47 + 0.20 


-1 


1341+0056 


18330 ± 436 


7.99 ± 0.10 


SP 


M7 


DAire 


18330 ± 608 


8.01 + 0.13 


-1 


1352+0910 


36154 ± 722 


7.49 ± 0.11 


SP 


M5 


DAire 


35331 ± 919 


7.31 + 0.16 


-1 


1415+0117 


32000± 2 ||0 


8.00 


p 


M6 


DA:ire 


73816 ± 4981 


8.43 + 0.18 


M0 


1448+0713 


12250 ± 256 


7.87 + 0.09 


SP 


M8 


DAire 


12110 + 264 


7.84 + 0.13 


M6 


1539+2706 


28000+ffg 


8.00 


p 


M9 


DA:ire 


36572 ± 447 


7.31 + 0.07 


-1 


1619+2533 


25595 ± 1139 


7.21 ± 0.15 


SP 


M5 


DAire 


25891 ± 1447 


7.30 + 0.20 


-1 


2255-0015 


17000+ffg 


9.oo±8;ll 


p 


Li: 


DA:ire: 


22550 ± 724 


7.88 + 0.12 


M6 


2331+1342 


11173 + 284 


8.79 + 0.15 


SP 


L6 


DAire 


18756 ± 1020 


7.96 ± 0.20 


-1 



